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I  Introduction 


This  technical  memorandum  discusses  the  performance  of  high  resolution 
methods  used  to  process  the  data  from  vertical  line  arrays  in  a  realistic  oceanic  sound 
speed  environment. 

Chapter  II  introduces  the  array  processing  structures  commonly  used  to  beam- 
form  line  arrays.  After  an  outline  of  the  frequency  wavenumber  representation  and  its 
associated  hypothesis,  the  conventional  beamformer  is  derived.  A  signal  model  that 
relates  the  covariance  matrix  to  the  power  and  the  directions  of  the  arrivals  is  introduced, 
and  the  beamforming  process  is  decomposed  into  two  steps.  First,  the  covariance  matrix 
is  estimated  and  a  measure,  based  on  the  covariance  matrix,  is  computed  to  yield  the 
amount  of  match  or  mismatch  between  the  field  and  a  replica  vector  that  characterizes  a 
particular  arrival  angle.  High  resolution  methods  are  introduced  and  a  well-known  adap¬ 
tive  beamforming  technique,  the  MVDR  processor,  is  described.  Some  performance  and 
implementation  issues  are  addressed.  Finally,  superresolution  methods  are  introduced 
and  the  MUSIC  eigenvector  method  is  presented. 

High  resolution  beamformers  are  well  known  to  fail  under  correlated  arrivals.  In 
Chapter  HI,  the  difference  in  the  structure  of  the  array  covariance  matrix  between  a 
scenario  with  correlated  arrivals  and  a  scenario  with  only  uncorrelated  arrivals  is 
presented.  The  observation  of  the  non-Toeplitz  character  of  the  covariance  matrix  for 
correlated  arrivals  leads  to  the  spatial  smoothing  techniques  used  to  decorrelate  corre¬ 
lated  arrivals.  These  preprocessing  techniques  are  described  and  the  reasons  how  and 
why  they  perform  are  discussed.  Then,  using  simulations  as  well  as  analytical  deriva¬ 
tions,  the  performance  of  modified  spatial  smoothing  is  compared  to  the  performance  of 
the  original  spatial  smoothing. 
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The  effects  of  wavefront  curvature  are  studied  in  Chapter  IV  in  the  context  of 
vertical  arrays  using  a  geometric  ray  theory  approach.  Based  on  Snell’s  law,  a  way  to 
generate  curved  wavefronts  is  derived  analytically.  The  importance  of  wavefront  curva¬ 
ture  is  first  evaluated  by  comparing  the  phase  between  curved  wavefronts  and  plane 
wavefronts,  both  corresponding  to  the  same  arrival  angles  for  a  realistic  scenario.  The 
behaviour  of  the  conventional,  MVDR,  and  MUSIC  beamformers  using  plane  wavefront 
replica  vectors  then  is  studied  in  a  curved  wavefront  environment.  The  largest  mismatch 
loss  occurs  for  near  horizontal  arrivals  and  the  aperture  length  is  not  critical  for  this 
rather  general  result.  Curved  wavefront  beamforming  then  is  proposed  and  the  behaviour 
of  the  processors  is  described  in  terms  of  bearing  responses  and  beam  patterns.  The  dis¬ 
tortion  caused  by  the  smoothing  transformation  is  studied  in  the  context  of  high  resolu¬ 
tion  methods.  An  approach  using  the  eigenvectors  of  the  smoothed  covariance  matrix  for 
each  replica  vector  allows  the  loss  due  to  mismatch  to  be  avoided.  Finally,  the  perfor¬ 
mance  of  the  MVDR  beamformer  after  spatial  smoothing  under  correlated  and  curved 
arrivals  is  discussed. 
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II  Array  Signal  Processing 


This  chapter  introduces  the  array  processing  structures  commonly  used  to  beam- 
form  line  arrays.  After  an  outline  of  the  frequency  wavenumber  representation  and  its 
associated  hypothesis,  the  conventional  beamformer  is  derived.  A  signal  model  that 
relates  the  covariance  matrix  to  the  power  and  the  directions  of  the  arrivals  is  introduced 
and  the  beamforming  process  is  decomposed  into  two  steps.  First,  the  covariance  matrix 
is  estimated  and  a  measure  based  on  the  covariance  matrix,  is  computed  to  yield  the 
amount  of  match  or  mismatch  between  the  field  and  a  replica  vector  that  characterizes  a 
particular  arrival  angle.  High  resolution  methods  are  introduced  and  a  well-known  adap¬ 
tive  beamforming  technique,  the  MVDR  processor,  is  described.  Some  performance  and 
implementation  issues  are  addressed.  Finally,  superresolution  methods  are  introduced 
and  the  MUSIC  eigenvector  method  is  presented.  These  three  processing  schemes  will  be 
used  in  the  subsequent  chapters. 


4 


n.l  Basic  Assumptions 

Array  processing  methods  decompose  the  pressure  field  into  a  directionality 
representation  by  assuming  homogeneous  wavefields  [Yenl977].  A  homogeneous 
wavefield  is  the  superposition  of  uncorrelated  elementary  waves  [Bohrnel987].  This 
mathematical  formulation  sets  a  clear  analogy  between  time  series  spectral  estimation 
and  array  processing  in  its  frequency-wavenumber  representation  of  the  random  field 
[Bohmel987,  Capon  1969] : 

P(t,?)  =  fexp[j  (o it  +  k’r)]dZ(a,k  )  (2.1.1) 

where  p  denotes  the  random  field  at  time  t  and  spatial  position  t ’,  w  the  circular  frequency 
in  rad/s  and  k  the  wavenumber  vector.  dZ  is  the  spectral  increment  of  the  random  spectral 
measure  Z(co ,  k)  with  covariance 

COv\dZ((o,k  )dZ  (a)',  k'  )1  =  8(a>-a>' )  5(  k-  k* )  S  (to  jT)  5a)  dk  (2.1.2) 

L  J  (2rt)4 

where  to  and  to'  are  circular  frequencies,  k  and  k’  wavenumber  vectors  ,  S  the  frequency- 
wavenumber  spectrum  of  the  random  field  p(t ,?)  and  8  the  dirac  delta  function. 

This  representation  is  a  generalization  of  the  Cramer  spectral  representation  for 
stationary  time  series  [Priestleyl981].  Its  main  assumption  is  that  spectral  increments  dZ 
at  two  distinct  frequencies  or  at  two  distinct  wavenumber  vectors  are  uncorrelated.  It 
constitutes  the  foundation  of  any  array  processing  method. 

Thus,  this  representation  is  one  where  the  receiving  array  measures  signals 
which  result  from  the  superposition  of  uncorrelated  elementary  waves.  These  elementary 
waves  generally  are  assumed  planar.  Such  a  homogeneous  wavefield  representation  is 
restrictive  since  it  excludes  correlated  arrivals  (e.g.  multipaths). 
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H.2  Conventional  Beamforming 


n.2.1  Lagged  Regression  Approach 

Conventional  beamforming  with  a  line  array  made  of  uniformely  spaced  sensors 
is  a  simple  illustration  of  the  already  mentioned  analogy  with  classic  time  series  spectral 
estimation  [Kayl988].  The  conventional  beamformer  is  the  simplest  processing  structure 
since  it  implements  a  linear  regression  on  the  sensor  outputs  either  in  the  time  domain  or 
in  the  frequency  domain  [Shumwayl988]. 

The  time  domain  implementation  is  known  as  delay-and-sum  beamforming.  A 
beam  is  formed  in  the  time  domain  by  summing  up  the  lagged  (delayed)  outputs  of  the 
array  sensors  [Johnson  1982,  De  Fattal988] 

y(/,0)=  Lamxm(t-xm)  (2.2.1) 

m=  0 

where  t  denotes  time,  x (t)  is  the  output  of  the  mth  sensor,  M  the  number  of  sensors, 
(am)m=o,M-i  the  shading  coefficients  and  (*„)«= o.m-i  the  delays  associated  with  each  sensor 
corresponding  to  a  particular  look  direction,  6.  The  time  delays  xn  are  determined  by 
making  additional  assumptions  on  the  propagating  signals.  If  one  considers  plane  waves 
impinging  on  a  line  array  with  equally  spaced  sensors,  the  delays  xm  are  given  by 


X-  = 


m  d  sin0 


(2.2.2) 


where  d  is  the  interelement  spacing,  c  the  sound  velocity  and  0  the  look  direction 
[DeFattal988].  The  delay-and-sum  beamformer  is  summarized  in  Figure  II.l. 

The  frequency  domain  regression  can  be  derived  by  taking  the  Fourier  transform 
of  the  output  y(t) 


Y(f,Q)=  l!  ane~J^  Xn(f)  (2.2.3) 

m=  0 

where  /  is  the  carrier  frequency,  Y(f,  0)  and  Xm(f)  the  Fourier  transforms  of  y(t)  and  x(t) 
respectively.  In  essence,  Equation  (2.2.3)  is  the  narrow-band  implementation  of  Equation 
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(2.2.1)  at  frequency  / .  Equation  (2.2.3)  can  also  be  expressed  as  an  inner  product : 


Y  (f  ,k  )  =  A"  X 

(2.2.4) 

[*.(/•)  *,(/•)  •••  X*-,(f)] 

(2.2.5) 

i7xfx-  a.c™" 

(2.2.6) 

the  superscript  H  denotes  the  Hermitian  operation  (complex  conjugate  transpose  opera¬ 
tion),  and  the  superscript  T  denotes  the  transpose  operation.  Assuming  plane  waves,  that 
is  using  Equation  (2.2.2)  and  letting 


k  =  *  “  (2.2.7) 

c 

Equation  (2.2.3)  is  in  the  form  of  a  spatially  windowed  Discrete  Fourier  Transform 
[Oppenheiml975]  : 


Y(fJc)=Zane-'2*kmXm(f)  (2.2.8) 

m=  0 

where  k  is  a  normalized  wavenumber. 

The  frequency-wavenumber  power  spectral  density  is  estimated  by  the  periodogram 

fp(f,k)= - - 'Y(fjc)\2  (2.2.9) 

La,2) 

where  M  is  the  number  of  sensors.  Ip  (f,k)  can  be  converted  into  a  power  spectral  den¬ 
sity  parametrized  by  the  arrival  angle  0  by  the  nonlinear  mapping  given  in  Equation 
(2.2.7)  [De  Fattal988].  The  narrow-band  conventional  beamformer  often  is  imple¬ 
mented  as  a  FFT  beamformer  [Williams  1968].  The  FFT  beamformer  is  summarized  in 
Figure  II.2.  Since  it  is  well  known  that  the  periodogram  gives  an  inconsistent  estimate  of 
the  power  spectral  density  [Priestley  1981],  additional  incoherent  averaging  over  succes¬ 
sive  time  snapshots  is  typically  performed.  The  estimate  then  becomes 

i  1  K-\  2 


P  ( /.*)= 


L/t'  *Ki‘° 


j  K- 1 

-F.Z  I  Y.ifjc)  i 


(2.2.10) 


5V) 

Af  i*o 
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where  YL(fJc )  is  FFT  beamformer  output  from  from  the  iA  time  snapshot. 


n.2.2  Formulation  in  Terms  of  the  Covariance  Matrix 

When  the  window  applied  is  the  rectangular  window,  the  a,  ’s  are  taken  equal  to 
unity  and  the  spectral  density  can  then  be  expressed  in  a  simplified  matrix  form  as  : 


iptf  ,k)  =  ±(EHX)(EHXf 


=  —ehxxhe 

M 


where  E  is  the  plane  wave  steering  vector. 


ET  _  j  e)  2ttk  ...  ej2xk(M-l) 


(2.2.11) 


(2.2.12) 


(2.2.13) 


and  X  XH  is  an  estimate  at  the  frequency  /  of  the  cross  spectral  density  matrix  based  on 
the  data,  or  covariance  matrix.  With  additional  averaging  over  time  snapshots,  the  con¬ 
ventional  beamformer,  also  called  the  Bartlett  beamformer,  is  given  by 


PB(f,k)=^-EH  RE 

M 


1 

R 

K  <  =  o 


(2.2.14) 


(2.2.15) 


where  X,  is  the  vector  of  Fourier  coefficients  corresponding  to  the  time  snapshot  /  .  R  is 
an  estimate  of  the  array  covariance  matrix.  When  a  weighting  window  is  appplied,  one 


W.*)=  T7  Ew  W  R  Ww  E 
Af 


where  the  (M,M)  diagonal  matrix  W  is 


yVT  =  ——=  [a0a,  aM., ] 


(2.2.16) 


(2.2.17) 


n  (element  index) 
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Figure  II. 2:  FFT  beamformer. 
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n.2.3  Matched  Filter  Interpretation 

Expressed  as  a  function  of  the  array  covariance  matrix,  the  conventional  or 
Bartlett  beamformer  allows  the  decomposition  of  the  processing  structure  into  two  com¬ 
ponents:  the  processor  itself  and  the  set  of  steering  vectors. 

Using  a  given  set  of  replica  vectors,  finding  directions  of  arrival  corresponds  to 
finding  the  angles  of  maximum  power  at  the  output  of  the  processor.  Therefore,  the 
beamforming  operation  can  be  interpreted  as  a  spatial  matched  filtering  operation.  Con¬ 
sidering  Equation  (2.2.12),  the  quantity  to  be  maximized  is 

Aw  XXH  A  =  I  Ah  X 1 2  (2.2.16) 

where  A  is  a  direction  vector  given  by  Equation  (2.2.6)  corresponding  to  an  arrival  angle 
0.  A  belongs  to  a  set  of  steering  vectors  or  replica  vectors  which  often  is  called  the  array 
manifold  [Schmidtl981,  Clarkel988].  The  choice  of  an  array  manifold  is  equivalent  to 
the  choice  of  a  wavefront  model.  Such  a  choice  is  guided  by  the  underlying  assumptions 
about  the  signal  propagation  and  wavefront  shape.  Generally,  the  plane  wave  array  man¬ 
ifold  is  selected. 

If  a  single  signal  impinges  on  the  array  at  a  given  but  unknown  angle  of  arrival, 
the  measured  array  signal  vector  X  given  by  Equation  (2.2.5)  will  be  an  amplitude 
weighted  plane  wavefront  vector.  From  the  Schwartz  inequality  the  output  of  the  beam- 
former  will  be  maximum  if  the  replica  vector  A  and  the  array  signal  vector  X  are  colinear. 
The  beamformer  provides  a  measure  of  the  mismatch  between  the  observed  field  and  a 
replica  vector  or  steering  vector.  A  match  implies  a  high  output  and  indicates  an  arrival 
direction  [Cox  1973]. 
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H.3  Covariance  Matrix  Representation 

In  the  previous  section,  it  was  shown  that  the  conventional  processor  operates  on 
the  array  covariance  matrix.  In  general,  any  beamformer  will  operate  on  the  array  covari¬ 
ance  matrix  since  it  contains  information  on  directionality.  Therefore,  it  is  of  interest  to 
discuss  the  covariance  matrix  representation. 

If  p  narrowband  signals  at  the  carrier  frequency  /  impinge  on  the  array,  the  time 
domain  array  signal  vector  at  the  input  of  the  beamformer  is  given  by  [Coxl973,  De 
Graafl985] 

r(0  =  Bs(f)  +  n(r)  (2.3.1) 

where  s (t)  is  the  (p,l)  signal  vector,  sometimes  called  the  envelope  vector  [Cad- 
zowl988],  n(r)  is  the  (Af  ,1)  vector  of  additive  noise  and  B  the  matrix  of  stacked  directions 
with  columns  equal  to 

B=[ao  A,  •••  A,_,]  (2.3.2) 

each  (M  ,1)  column  vector  A ,  corresponds  to  an  arrival  angle  0,.  The  equivalent  narrow- 
band  frequency  domain  representation  is,  at  the  signal  carrier  frequency  / , 

r  =  Bs  +  n  (2.3.3) 

where  s  is  the  (p,l)  signal  vector,  n  is  the  additive  noise  (Af  ,1)  vector,  and  B  the  matrix  of 
stacked  directions.  The  corresponding  model  covariance  matrix  Rmad,1  =  r  r11  is  thus 

Rmo*,=BSBw+Q  (2.3.4) 

where  S  is  the  (p,p)  signal  covariance  matrix  given  by  £(ssw)  and  Q  is  the  noise  covari¬ 
ance  matrix  given  by  £(n  nH)  with  £  denoting  the  expectation. 

Under  the  assumption  of  a  homogeneous  wavefield,  the  signals  are  uncorrelated, 
and  one  can  show  that  the  (p,p)  signal  matrix  S  is  of  full  rankp.  Recalling  that  S  =  £(ssw) 
with  sr  =  [j0*i  •••  5p_i],  uncorrelated  arrivals  means  that  the  signal  envelopes  are  such 
that  for  i*j. 
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£  fa  sj)  =  0  (2.3.5) 

where  *  denotes  the  complex  conjugate  operation.  If  s,  and  sj  are  correlated,  then 

E(Sis‘)*0  (2.3.6) 

Restricting  the  discussion  to  two  arrivals,  the  signal  covariance  matrix  is 


f£(i0iS)  £(*0*1)1 

S  =  •  .  (2.3.7) 

E(s0  si)  £(*i  si) 

k  * 

with  E(sjs’)  =  a ?,  the  power  of  the  iA  signal,  and  for  i  *  j,  E(s,  s‘)  =  aiajplj,  where  p,^  =  p*, 
is  the  complex  correlation  given  by 


_  E  & 

PiJ~  £(l*,l2)£(l*>l2) 

Then,  the  signal  covariance  matrix  S  can  be  written  [Reddy  1987] 

a02  o0  o,p 

S  —  •  j 

<T0OiP  <Tf 


(2.3.8) 


(2.3.9) 


For  fully  coherent  arrivals  I  p  I  =  1,  hence  det  S  =  ai  a,2  (1  -  t  p  1 2)  =  0. 

The  standard  model  for  the  noise  covariance  Q  is  a  sensor  noise  uncorrelated 
from  sensor  to  sensor  [Nickell988,  Schmidtl986].  Under  this  assumption,  the  noise 
covariance  matrix  is  Q  =  a2!,  where  a2  is  the  noise  power  and  I  the  identity  matrix.  It  fol¬ 
lows  that  the  (M ,  M )  covariance  matrix  is  of  full  rank. 

Furthermore,  if  one  assumes  the  elementary  waves  of  the  field  to  be  planar,  the 
covariance  matrix  which  is  Hermitian  by  construction  (i.e.  R H  =  Rwii  ),  also  is 

Toeplitz  (constant  along  the  diagonals). 

Rmo&z  =  ^o2  Ao  A0W  +  CTj2  Aj  AjW  +  +  Op-2  Ap_i  Ap_]W  (2.3.10) 

Since  A,  =  £  l  eJ2xK  e'2*2*  ■  ■  •  eJ2xtM~1)'  j ,  each  outer  product  in  the  sum  given  by  Equation 
(2.3.9)  is  Toeplitz : 


A,  A,  "  = 


1  e-'2* 

...  e-/2*<M-l)fc- 

e*7*  1 

• 

1 

eJ2K<M~\)k  g J 2ic<M -2)fc 

1 

(2.3.11) 
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so  that  R mm  is  both  Toeplitz  and  Hermitian.  The  Toeplitz  character  of  the  array  covari¬ 
ance  matrix  is  a  result  of  the  duality  between  homogeneous  plane  wavefields  and  wide 
sense  stationary  time  series. 


H.4  Adaptive  Beamforming 

n.4.1  Derivation  of  the  MVDR  beamformer 

The  narrow-band  beamforming  problem  can  be  viewed  as  finding  a  set  of  com¬ 
plex  weights  or,  equivalently,  a  steering  vector  given  by  Equation  (2.2.6)  [De  Graafl985] 
corresponding  to  the  signal  direction  of  arrival.  The  plane  wave  assumption  only 
corresponds  to  a  particular  signal  model.  The  conventional  technique  is  the  most  widely 
used,  the  best  known  and  the  most  robust,  but  it  has  the  least  resolution.  In  an  effort  to 
improve  the  resolution,  data  adaptive  techniques  can  be  used.  The  increase  in  resolution 
results  from  the  full  use  of  the  information  summarized  in  the  covariance  matrix,  along 
with  specific  assumptions  about  the  signal  model 

An  adaptive  beamformer  determines  a  data  dependent  set  of  weights  as  the  solu¬ 
tion  of  a  constrained  optimization  problem  [Griffiths 1977,  Monzingol985,  Van 
Veenl988].  Thus,  it  is  referred  to  as  an  optimum  beamformer  in  the  literature  [Coxl973, 
Reddy  1987]. 

The  classical  optimization  criterion  is  to  minimize  the  output  power  or  variance 
under  the  constraint  of  passing  the  look  direction  signal  undistorted  [Coxl973, 
Kanasewichl975,  Marplel987,  Lacossl971] : 

MIN  [ww  R  w]A.w=]  (2.4.1) 

where  W  is  the  unknown  weight  vector,  A  the  look  of  direction  vector  or  replica  vector 
from  the  array  manifold,  as  in  Equation  (2.2.6),  and  R  the  cross- spectral  density  matrix 
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estimated  from  the  data  given  by  Equation  (2.2.5).  As  before,  one  asssumes  that  the  data 
have  been  Fourier  analysed  and  one  considers  a  single  frequency  /.  Using  the 
Lagrange’s  multiplier  method  [Lacossl971,  Kanasewitchl975,  Marplel987],  the  weight 
vector  W  is  given  by 

W  =  — (2.4.2) 
A H  R->  a  K  ’ 

and  the  wavenumber  spectrum  is 

P  uvdr  (f » k )  =  H  ,  (2.4.3) 

where  M  is  the  number  of  sensors.  This  is  the  so-called  Maximum  Likelihood  method  of 
Capon  [Caponl969],  which  is  the  maximum  likelihood  estimator  in  the  case  of  Gaussian 
noise  (i.e.  it  minimizes  the  log- likelihood  under  the  assumption  of  Gaussian  statistics). 
This  beamformer  also  is  called  the  Minimum  Energy  [De  Graafl985]  or  Minimum  Vari¬ 
ance  Distortionless  Response  (MVDR)  beamformer  [Zoltowskil988].  The  MVDR 
beamformer  stands  as  a  high  resolution  method  compared  to  the  conventional  beamform¬ 
ing.  For  this  reason,  it  is  used  in  various  areas  such  as  the  seismic  field  [Capon  1969,  Bag- 
geroerl974,  Baggeroerl982,  Hsu  1986],  the  oceanographic  field  [01tman-shayl984]  as 
well  as  the  underwater  acoustics  arena  [Owsleyl985]. 

n.4.2  Performance  of  the  MVDR  Beamformer 

The  performance  of  the  MVDR  beamformer  has  been  studied  rather  extensively 
[Coxl973,  De  Graafl985].  The  main  results  are  presented  here  succintly.  In  the  context 
of  matched  filtering,  the  MVDR  processor  is  such  that  a  slight  mismatch  in  the  replica 
vector  results  in  serious  signal  reduction  at  its  output  [Cox  1973].  For  the  adaptive  array 
pointing  in  a  given  direction,  the  signal  arriving  in  another  bearing  is  an  interference  to 
be  cancelled.  The  stronger  the  interference,  the  more  suppression  the  designed  filter  is 
able  to  put  on  it.  This  suppression  is  based  on  the  covariance  structure  across  the  array 


that  samples  the  field.  Depending  on  the  signal-to-noise  ratio  SNR ,  the  resolution  varies. 
For  a  line  array  and  two  signals  30  dB  above  the  noise  level,  the  resolution  of  the  MVDR 
beamformer  [Coxl973]  is  slightly  over  three  times  the  Rayleigh  resolution  of  conven¬ 
tional  processing  [Burdicl984]  which  is  given  in  radians  by 

e  =  j  (2.4.4) 

where  X  is  the  acoustic  wavelength  and  L  is  the  length  of  the  array.  When  SNR  -» 0,  the 
resolution  of  the  MVDR  processor  decreases  and  becomes  comparable  to  the  conven¬ 
tional  processor  one  [Coxl973,  Kanasewichl975,  Zoltowskil988].  This  improved  reso¬ 
lution  results  in  a  higher  sensitivity  to  mismatch.  Some  cost  may  thus  be  associated  to 
this  higher  sensitivity  to  mismatch  in  terms  of  signal  to  noise  ratio  improvement  or  pro¬ 
cessing  gain.  If  the  array  manifold  is  chosen  with  improper  assumptions  about  the  true 
underlying  wavefronts,  the  resolving  capability  of  the  MVDR  beamformer  can  be  less 
than  that  of  the  conventional  processor  [Seglisonl970,  McDonough 1972]. 

n.4.3  Implementation  Issues 

The  MVDR  beamformer  is  implemented  by 

(1)  performing  the  Fourier  analysis  at  the  output  of  each  sensor  using  conventional 
FFT  based  techniques  on  time  subsegments 

(2)  computing  the  cross- spectral  matrix  based  on  the  data,  using  Equation  (2.2.5) 

(3)  inverting  the  cross-spectral  matrix 

(4)  computing  the  angular  wavenumber  spectrum  given  by  Equation  (2.4.3). 

It  is  assumed  that  the  signals  are  superimposed  on  an  incoherent  noise.  A  nonne¬ 
gative  definite  cross-spectral  matrix  is  ensured  theoretically  by  averaging  at  least  M  dyad 
products  X,  X,w,  each  dyad  product  corresponds  to  a  dimension.  Averaging  K  cross- 
spectral  matrices  as  in  Equation  (2.2.15),  corresponds  to  averaging  across  K  time  seg- 
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ments,  each  X;  results  from  a  time  subsegment.  Even  when  K  >M,  the  matrix  inversion 
requires  special  attention,  that  is  double  precision  arithmetic  and  an  additional  check  for 
matrix  ill-conditioning.  A  cross-spectral  matrix  based  on  the  data  often  is  numerically 
noninvertible  ;  this  situation  is  characterized  by  the  ratio  of  its  smallest  to  its  largest 
eigenvalue  with  a  practical  criterion  being  a  ratio  below  10-6.  The  solution  [Caponl969] 
then  is  to  add  a  small  amount  of  incoherent  noise  to  the  cross-spectral  matrix  R,  that  is  to 
add  to  the  main  diagonal  of  R  : 


e  =  y 


tr( R) 
M 


(2.4.5) 


where  tr  denotes  the  trace  operation,  and  y  the  fraction  of  noise,  and 


represents  the 
M 


mean  power  across  the  sensor  array.  Typical  values  for  y  are  in  the  range  of  10"*  to  10'2. 
Taking  y=  KT4  corresponds  to  introducing  an  uncorrelated  sensor  noise  40  dB  below  the 
average  sensor  power  level. 

This  stabilization  procedure  is  equivalent  to  increasing  the  incoherent  sensor 
noise  and  thus  produces  a  small  bias  in  the  power  spectral  density  estimate,  in  addition  to 
the  one  due  to  the  data  sequence  length  limitation  which  is  given  by  [Capon  1970] 


bias (dB)=  10 logK  +  1  (2.4.6) 

K 

The  asymptotic  distributions  (i.e.  for  very  long  data  sequences)  of  the  MVDR  spectrum, 
pmvdr  tf.*),  is  a  Wishart  distribution,  or  a  multidimensional  chi-square  [Caponl970]  with 
a  number  of  degrees  of  freedom  given  by 


d.o.f.  =  2(K-M +1) 


(2.4.7) 
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H.5  Eigenvector  Beamforming 

The  direction  of  arrival  estimation  problem  can  be  solved  using  spectral  estima¬ 
tion  techniques.  However,  sufficient  angular  resolution  may  not  be  achieved  with  con¬ 
ventional  or  adaptive  techniques.  Direction  of  arrival  or  frequency  estimators  have  been 
developed  to  yield  maximum  resolution  based  on  the  covariance  matrix  eigenstructure. 
These  alternative  methods  provide  super  resolution  but  not  an  estimation  of  the  power 
spectral  density. 

The  eigenvector  beamformer  has  the  same  assumptions  as  the  adaptive  beam- 
fonner,  the  covariance  matrix  is  the  superposition  of  a  signal  covariance  matrix  and  a 
noise  covariance  matrix,  as  in  Equation  (2.3.4).  There  are  p  signals  impinging  on  the 
array  and  the  noise  is  a  sensor  noise  uncorrelated  from  sensor  to  sensor.  The  p  largest 
eigenvalues  of  the  covariance  matrix  allow  the  separation  of  the  p  signals  from  the  noise, 
the  corresponding  eigenvectors  span  the  signal  subspace.  The  other  M  -  p  eigenvectors 
corresponding  to  the  lowest  eigenvalues  span  the  noise  subspace  which  is  orthogonal  to 
the  signal  subspace  [Schmidtl981].  This  geometrical  formulation  leads  to  two  classes  of 
estimators,  one  based  on  the  signal  subspace  and  the  other  based  on  the  noise  subspace 
[Kayl988,  Marplel987]. 

The  MUSIC  (Multiple  Signal  Classification)  algorithm  [Clarkel988,  Can- 
tonil980]  is  a  principal  candidate  among  the  various  spectral  eigenvector  estimators,  and 
is  widely  referenced  in  the  literature.  MUSIC  is  a  noise  subspace  estimator  where  the 
orthogonality  of  the  noise  and  signal  subspaces  plays  a  major  role.  This  estimator  is 
given  by 


(2.5.1) 


ruvs,c  ~ 

where  A  is  the  steering  vector  as  in  Equation  (2.1.5)  and  R*  the  noise  only  covariance 
matrix,  given  by 
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U-p- 1 

Rjv=  I  \kVkH  (2.5.2) 

k=C 

where  Vt  is  the  noise  eigenvector  of  the  data  covariance  matrix  given  by  Equation 
(2.2.15).  V*  is  such  that  V*w  V*  =  l. 

MUSIC  is  implemented  by 

(1)  the  eigenvalue  and  eigenvector  decomposition  of  the  covariance  matrix 

(2)  the  selection  of  the  number  of  signals  p 

(3)  the  computation  of  the  noise  only  covariance  matrix  given  by  Equation  (2.5.2) 

(4)  the  computation  of  the  MUSIC  spectrum  given  by  Equation  (2.5.1)  using  the 
plane  wave  array  manifold 

The  MUSIC  spectrum  is  the  inverse  of  the  squared  distance  of  the  projection  of 
A  on  the  noise  subspace.  When  scanning  all  the  possible  replica  vectors  of  the  array  man¬ 
ifold,  global  minima  occur  if  this  distance  is  zero,  or  equivalently,  when  A  belongs  to  the 
signal  subspace.  Then,  the  MUSIC  spectrum  theoretically  will  show  a  peak  to  infinity.  In 
practice,  the  peaks  are  finite  because  one  has  only  an  estimate  of  the  covariance  matrix. 

The  noise  only  covariance  matrix  R*  is  of  rank  M  -p,  thus  rank  deficient.  The 
quadratic  form  \hrn\  in  Equation  (2.5.1)  is  positive  semi-definite.  The  noise  covariance 
matrix  has  M  -p  eigenvalues  equal  or  close  to  unity,  and  p  eigenvalues  close  to  zero.  To 
avoid  numerical  problems  when  computing  the  quadratic  form  in  Equation  (2.5.1),  (i.e., 
small  negative  values),  a  small  fraction  of  noise  can  be  added  to  the  main  diagonal  of  RN , 
typical  fractions  of  noise  are  between  10"8  and  lCr6.  A  fraction  of  noise  equal  to  10~7 
corresponds  to  clipping  the  direction  of  arrival  peaks  to  approximately  70  dB. 

The  angular  resolution  enhancement  of  the  eigenanalysis  based  methods 
[Kayl988,  Marplel987]  and  more  precisely  of  MUSIC  [Schmidt86,  Nickell987]  is 
recognized  widely.  An  explicit  relationship  between  the  MVDR  beamformer  and  the 
MUSIC  algorithm  is  derived  in  [Nickell988]  where  the  MVDR  beamformer  is  shown  to 
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be  a  weighted  projection  on  the  noise  subspace,  like  MUSIC,  when  the  signal  to  noise 
ratio  goes  to  infinity.  This  explains  the  superior  resolution  of  MUSIC  with  respect  to  the 
MVDR  beamformer. 

One  of  the  major  tasks  in  practice  is  to  determine  the  number  of  signals  imping¬ 
ing  on  the  array  [Schmidtl981,  Waxl984].  Three  methods  can  be  used 

(1)  the  number  of  signals  is  determined  by  observing  the  spread  of  the  eigenvalues 
of  the  covariance  matrix  and  using  a  threshold  level  to  seperate  high/signals  from 
low/noise  eigenvalues  [Waxl984,  Oswleyl988] 

(2)  the  number  of  signals  is  determined  using  hypothesis  testing  on  the  multiplicity 
of  the  smallest  eigenvalues  [Schmidtl980,  Waxl984] 

(3)  information  theoretic  criteria  of  model  selection  allow  the  determination  of  the 
number  of  signals  without  any  subjective  judgement  as  in  (1)  and  (2)  [Waxl985] 
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HI.  Processing  Correlated  Arrivals 

High  resolution  beamformers  are  well  known  to  fail  under  correlated  arrivals 
[Whitel979,  Cantonil980].  Here,  the  difference  in  the  structure  of  the  array  covariance 
matrix  between  a  scenario  with  correlated  arrivals  and  a  scenario  with  only  uncorrelated 
arrivals  is  presented.  The  observation  of  the  non-Toeplitz  character  of  the  covariance 
matrix  for  correlated  arrivals  leads  to  the  smoothing  techniques  [Shanl985,  Willi- 
amsl988]  used  to  decorrelate  correlated  arrivals.  These  preprocessing  techniques  are 
described  and  the  reasons  how  and  why  they  perform  are  discussed.  Then,  the  perfor¬ 
mance  of  modified  spatial  smoothing  is  compared  to  the  performance  of  the  original  spa¬ 
tial  smoothing  [Reddyl987],  using  simulations  as  well  as  analytical  derivations. 
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in.l  Limitation  of  the  High  Resolution  Methods 

The  analogy  between  time  series  analysis  and  spatial  array  processing  leads  to 
the  formulation  of  directionality  in  the  framework  of  homogeneous  wavefields. 
Nevertheless,  array  processing  is  not  identical  to  time  series  power  spectral  density  esti¬ 
mation.  An  oceanic  environment  is  bounded  at  the  surface  and  the  bottom  and  these 
boundaries  are  known  to  have  a  major  impact  on  the  pressure  field  [Yenl977],  Multipath 
propagation  produces  correlated  arrivals  at  the  receiving  line  array  [Urickl983]. 

Techniques  such  as  the  adaptive  and  eigenvector  methods  are  high  resolution 
and  lack  robustness  as  soon  as  the  structure  of  the  array  covariance  matrix  significantly 
departs  from  the  assumed  signal  model.  Under  such  conditions,  they  generally  report 
improper  arrival  directions  or  incorrect  power  levels. 

The  severe  limitation  of  the  performance  of  the  MVDR  beamformer  under  corre¬ 
lated  arrivals  are  well  known  [Whitel979,  Cantonil980].  The  adaptive  filter  uses  the 
non-look  direction  signals  to  minimize  its  output  power  and  since  the  signals  are  corre¬ 
lated  with  the  look  direction  signals,  the  latter  is  suppressed.  This  process  is  the  so-called 
signal  cancellation  phenomena  in  adaptive  antennas  [Widrowl982].  On  an  algorithmic 
level,  the  adaptive  technique  fails  because  the  model  signal  covariance  is  singular 
[Gabriell980]. 

The  MUSIC  algorithm  also  does  suffer  from  arrival  intercorrelation 
[Schmidtl981].  If  there  are  two  correlated  arrivals  on  the  array,  the  signal  subspace  is 
degenerate  of  rank  one:  the  direction  of  arrivals  are  associated  to  the  same  single 
waveform  corresponding  to  a  single  emitter  [Gabriel  1980]. 

From  a  practical  standpoint,  the  high  resolution  MVDR  and  MUSIC  methods  are 
unable  to  resolve  correlated  arrivals  as  indicated  by  the  results  in  [White  1979]  and  [Wil- 
liamsl988].  This  point  is  further  illustrated  by  Figure  III.  1  and  IEI.2. 
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Figure  m.l:  MVDR  beamforming  on  3  plane  wave  arrivals  at  -40°,  -35"  and  -30°. 
Panel  A:  the  arrivals  are  uncorrelated,  Panel  B:  the  arrivals  at  -  V  and  -30°  are  fully 
correlated  (p  =  1).  The  array  has  32  sensors  with  25  m  spacing  and  operates  at  20  Hz  in 
an  isovelocity  medium  of  sound  speed  1498  m/s. 
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Figure  m.2:  MUSIC  beamforming  on  3  plane  wave  arrivals  at  -40°,  -35®  and  -30". 
Panel  A:  the  arrivals  are  uncorrelated,  Panel  B:  the  arrivals  at  -35"  and  -  30"  are  fully 
correlated  (p  =  1).  The  array  has  32  sensors  with  25  m  spacing  and  operates  at  20  Hz  in 
an  isovelocity  medium  of  sound  speed  1498  m/s. 
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III.2  Structure  of  the  Covariance  Matrix  Under  Correlated  Arrivals 

It  was  shown  earlier  that  full  coherence  among  arrivals  is  equivalent  to  a  rank 
deficiency  in  the  signal  covariance  matrix.  This  explains  why  the  adaptive  and  eigenvec¬ 
tor  methods  fail  under  correlated  arrivals.  Since  the  available  data  is  the  array  covariance 
matrix,  the  structure  of  the  covariance  matrix  for  an  inhomogeneous  field  is  now  investi¬ 
gated. 

It  was  shown  earlier  that  the  array  covariance  matrix  of  a  homogeneous  field 
using  the  plane  wave  array  manifold  is  both  Toeplitz  and  Hermitian.  In  a  realistic 
environment,  arrivals  such  as  multipaths  are  correlated  so  that  the  field  is  not  homogene¬ 
ous  (i.e.  not  stationnary  in  space).  Thus,  although  it  is  Hermitian  by  construction,  the 
covariance  matrix  of  a  line  array  is  not  Toeplitz  with  the  plane  wave  array  manifold 
[Gabriel  1980].  The  model  covariance  matrix  for  two  correlated  arrivals  can  be  expressed 
as  [Reddy  1987] 

Oo  CToOipl  rAw| 

R*od*/ =  [Ao  Aj]  •  2  ah  +<^I  (3.2.1) 

<J0OiP  Oi 

where  A,  is  the  direction  of  the  i1*  arrival,  a?  the  power  of  the  iA  arrival,  and  p  the  com¬ 
plex  correlation  between  the  two  arrivals,  I  p I  si,  and  o2  the  uncorrelated  sensor  noise 
power.  Using  the  plane  wave  array  manifold,  one  has  with  a  M  element  array 

A,=[l^  e*™-'*]7  (3.2.2) 

where  is  the  normalized  wavenumber  of  the  i*  arrival. 

The  upper  diagonal  terms  are  given  by 

Rmo*i  O'  i  +  1)  =  <To  e-'2**  +  a?  e~i2xk'  + 

o0a,|[p  *>»«*•-*•>)  e~j2xk'+  [p*  *>»<*« e~J2Kk'  ■  (3.2.3) 

for  i  =  l.Af-l.  Inspecting  these  terms,  one  concludes  that  the  matrix  is  not  Toeplitz.  The 
last  term  due  to  correlation  (p*0)  is  responsible  for  that  character.  Correlation  can  be 
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thought  of  as  introducing  a  modulation  term  in  the  diagonals  of  the  array  covariance 
matrix. 

The  effects  of  correlation  on  the  covariance  matrix  can  be  observed  on  Figure 
m.3  and  Figure  EII.4  by  comparing  the  magnitude  and  phase  of  the  upper  diagonals  of 
the  array  covariance  matrix  for  a  field  corresponding  to  two  uncorrelated  and  two  fully 
correlated  arrivals  (p=  l).  The  order  of  the  diagonals  can  be  identified  by  their  length. 

The  arrivals  are  unit  power  plane  waves  impinging  on  a  10  j  wavelength  line  array  of  32 

sensors,  at  -  10°  and  -  20°  incidence  angles.  The  interelement  spacing  is  25  m,  the  operat¬ 
ing  frequency  is  20  Hz  and  the  sound  speed  is  1498  m/s.  The  magnitude  and  phase  in  the 
uncorrelated  case  of  Figure  M.3  are  constant  on  each  diagonal,  while  they  both 
significantly  vary  in  the  correlated  case  on  Figure  M.4. 

This  observation  of  this  non-Toeplitz  character  is  important  since  it  points  out 
the  difference  in  nature  between  uncorrelated  arrivals  and  correlated  arrivals,  and  leads  to 
the  processing  methods  designed  to  limit  the  negative  effects  of  correlation. 
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Figure  m.3:  Magnitude  and  phase  of  the  upper  diagonals  of  the  32  by  32  covariance 
matrix.  Diagonals  are  identified  by  looking  at  their  length.  Matrix  index  corresponds  to 
element  index  along  the  diagonal.  The  32  sensor  line  array  (25  m  interelement  spacing) 
receives  two  uncorrelated  arrivals  at  - 10®  and  -  20"  incidence  angles  at  20  Hz.  The  sound 
speed  is  1498  m/s. 
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Figure  m.4:  Magnitude  and  phase  of  the  upper  diagonals  of  the  32  by  32  covariance 
matrix.  Diagonals  are  identified  by  looking  at  their  length.  Matrix  index  corresponds  to 
element  index  along  the  diagonal.  The  32  sensor  line  array  (25  m  interelement  spacing) 
receives  two  correlated  arrivals  (p  =  1)  at  -  10*  and  -  20*  incidence  angles  at  20  Hz.  The 
sound  speed  is  1498  m/s. 
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m.3  Processing  Methods 

The  methods  which  limit  the  effects  of  spatial  correlation  operate  on  the  covari¬ 
ance  matrix  to  make  it  "more  Toeplitz"  by  performing  some  kind  of  averaging  [Line- 
bargerl988].  The  straight  arithmetic  and  geometric  averaging  of  the  covariance  matrix 
diagonals,  based  on  the  concept  of  redudancy  [Linebargerl988],  have  been  proposed  in 
the  literature  [Naidul988,  Hsul986].  However,  forcing  a  Toeplitz  structure  does  not  sys¬ 
tematically  result  in  a  non  negative  definite  cross-spectral  matrix  and  is  considered  not 
viable  [Linebargerl988].  The  major  averaging  techniques  proposed  in  the  literature  are 

(1)  Frequency  averaging  which  is  averaging  covariance  matrices  at  different  fre¬ 
quencies  [Wangl985,  Zhul988].  This  method  is  relevant  only  to  broadband 
situations. 

(2)  Spatial  Averaging  or  Spatial  Smoothing,  developed  for  narrowband  problems 
(and  for  a  plane  wavefront  array  manifold)  [Gabriell980,  Shanl985,  Sul986, 
Reddy  1987,  Takaol987,  Linebargerl988].  This  method  is  particularly  relevant 
to  the  problems  discussed  here. 

III.4  Spatial  Smoothing 

III.4.1  Description  and  Interpretation 

The  concept  of  spatial  smoothing  comes  form  the  idea  that  the  relative  phase 
between  two  correlated  arrivals  changes  in  space  [Gabriell980,  Widrowl982],  In  prac¬ 
tice,  moving  the  array  spatial  position  is  out  of  question.  Spatial  smoothing  implements 
this  idea  by  averaging  covariance  matrices  estimated  on  subarrays  extracted  from  the  full 
aperture  and  corresponding  to  slightly  different  spatial  position.  This  averaging  induces 
a  reduction  of  the  terms  due  to  correlation  [Linebargerl988]  because  the  phase  variations 
are  eventually  averaged  to  zero. 

The  technique  uses  the  special  geometry  of  line  arrays  with  equi-spaced  sensors 
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[Linebargerl988]  and  computes  the  average  of  the  covariance  matrix  estimated  on  subar¬ 
rays,  where  each  subarray  shares  all  but  one  of  its  sensors  with  an  adjacent  subarray 
[Shanl985].  If  the  full  aperture  array  has  M  sensors,  and  the  subanrays  have  5  sensors, 
this  is  equivalent  to  averaging  M  -  s  + 1  blocks  of  s  by  s  element  extracted  from  the  full 
array  M  by  M  covariance  matrix  along  the  main  diagonal  [Takaol987].  This  procedure  is 
summarized  in  Figure  III. 5.  It  can  also  be  viewed  as  a  low-pass  filtering  operation  [Line- 
bargerl988].  Furthermore  averaging  decreases  the  variance  of  the  estimate  of  the  cross- 
spectral  matrix. 

Spatial  smoothing  is  implemented  as  follows  [Shanl985].  If  X  if  the  full  aper¬ 
ture  (A/,1)  array  vector,  and  X,  the  iA  (5,1)  subaperture  signal  vector 

X,T  =  [x,(f )  XM(f)  ■  ■  ■  X,«.(o]  (3-4.1) 

the  spatially  smoothed  covariance  matrix  is 

Ms 

=  T7-7-T  E  x‘  X-W  (3-4-2) 

Af  —  5  +  1 1  =  o 


ni.4.2  Effects  of  Spatial  Smoothing 

It  is  shown  in  [Shan  1985]  that  if  p  correlated  signals  impinge  on  the  array,  at 
least  p  averages  over  subarrays  of  p  elements  must  be  performed  in  order  to  recover  the 
rank  of  the  signal  covariance  matrix,  that  is  to  have  det  S  *  0  (i.e.  S  is  a  p  by  p  matrix  of 
rank  p).  This  result  is  of  great  importance  because  it  indicates  the  number  of  averages 
necessary  to  allow  the  determination  of  the  direction  of  arrivals  with  high  resolution 
methods  such  as  the  MVDR  beamformer  or  the  MUSIC  algorithm.  It  leads  to  the  rule  of 

thumb  for  spatial  smoothing  that  a  M  sensor  array  is  able  to  detect  at  most  -y  correlated 

arrivals  [Shanl985]. 

In  addition  to  their  direction  of  arrival,  it  often  is  of  interest  to  estimate  the 
power  of  the  arrivals.  The  rule  of  thumb  described  earlier  is  of  little  value  with  respect  to 
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this  problem  since  a  full  rank  matrix  does  not  ensure  total  decorrelation  of  the  arrivals, 
that  is  a  diagonal  signal  covariance  matrix.  Partial  correlation  among  the  arrivals  leads  to 
partial  signal  cancellation  for  the  adaptive  beamformers. 

The  spatial  smoothing  decorrelation  rate  for  two  correlated  arrivals  is  studied  in 
[Reddyl987].  It  is  shown  that  the  model  covariance  matrix  corresponding  to  the  ikrt 
subarray  can  be  expressed  as 


R*  =  B  S  (&~lf  Bh  +  a2 1  (3.4.3) 

where  R*  is  the  subarray  covariance  matrix  and  corresponds  to  Xt  x",  B  is  the  stacked 
plane  wave  direction  matrix  (a  s  by  d  matrix,  where  d  is  the  number  of  signals),  S  is  the  d 
by  d  signal  covariance  matrix  and  «t>  the  diagonal  matrix  given  by 
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(3.4.4) 


where  x,  represents  the  elementary  time  delay  between  two  consecutive  sensors  for  the  i,h 
arrival. 

This  formulation  shows  that  the  effects  of  spatial  smoothing  can  be  interpreted  in 
a  useful  way  in  terms  of  a  modified  signal  covariance  matrix,  the  smoothed  covariance 
matrix: 


S=  -p-  I  S ( Y* 

A  k  —  1 


(3.4.5) 


where  K  is  the  number  of  averages,  also  called  smoothing  index.  Restricting  the  discus¬ 
sion  to  two  correlated  arrivals,  [Reddy  1987]  indicates  that  perfect  decorrelation  is 
achieved  when  the  ( 1 ,2)  term  of  S  is  equal  to  zero.  This  term  is  given  by 


-  GnCTiP  K 

S(l,  2)  =  — p—  I  exp 

A  k  =  1 


-j(k- l)2rc—  (sin 0O -  sinGj) 


(3.4.6) 
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S(l.  2)  = 


CT1CT2P 


1  -  exp 


2nKA  ,  .  _  .  a  . 

— - — (sin80  -  sin0,) 


K 


1  -  exp 


■^-(sin90-sin8i) 


(3.4.7) 


where  p  is  complex  correlation  between  the  two  arrivals  of  power  Co,  a?  and  physical 
arrival  angles  80  and  0lt  A  is  the  array  interelement  spacing,  A.  the  acoustic  wavelength. 

The  minimum  value  of  the  smoothing  index  necessary  to  completely  decorrelate 
the  two  arrivals  is 


K  = - - -  (3.4.8) 

A  I  sin  0O  -sin0, 1  v  ' 

The  values  of  the  smoothing  index  K  are  tabulated,  for  an  array  with  half  wavelength  ele¬ 
ment  spacing,  as  a  function  of  the  angular  pair  separation  1 80  -  0i  l  and  the  position  of  the 
pair  referenced  by  1 80 1  (with  1 80 1  <  1 81 1 )  in  Table  IE.  1 .  An  arrival  pair  with  5°  separation 
arriving  at  10°  requires  at  least  12  averages  to  achieve  total  decorrelation  or,  equivalently, 
allow  the  MVDR  beamformer  to  yield  the  proper  powers. 

Equation  (3.4.8)  shows  that  the  number  of  averages  necessary  to  decorrelate  two 
close  arrivals  greatly  depends  on  the  pair  angular  spacing.  Two  closely  spaced  arrival 
are  very  difficult  to  decorrelate  [Reddyl987].  In  practice,  full  decorrelation  of  the 
arrivals  may  not  be  possible  if  the  array  has  a  limited  number  of  sensors.  The  effective 
aperture  is  reduced  from  M  to  s  which  results  in  lower  resolution.  Although  decorrela¬ 
tion  is  achieved  at  the  cost  of  resolution  and  numerical  processing  burden,  spatial 
smoothing  is  necessary  to  use  with  high  resolution  beamformers. 
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Figure  in.5:  Summary  of  spatial  smoothing  [Takaol987]. 
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Table  III.  1 :  Smoothing  index  as  a  function  of  arrival  pair  spacing  (lines  in  Table)  and 
position  from  broadside  (0°).  The  array  has  a  lA  wavelength  interelement  spacing. 
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IIL4.3  Simulation  Results 

m.4.3.1  Effects  on  the  Array  Covariance  Matrix 

Spatial  smoothing  is  a  preprocessing  technique  that  operates  on  the  full  array 
covariance  matrix  to  obtain  a  "more  Toeplitz"  matrix.  This  preprocessing  controls  the 
performance  of  the  high  resolution  techniques  by  its  ability  to  decorrelate  the  arrivals. 

To  illustrate  the  process  of  decorrelation  achieved  by  spatial  smoothing,  the 
diagonals  of  the  array  covariance  of  a  field  composed  of  two  correlated  arrivals  are 
inspected  under  several  smoothing  scenarios.  As  in  Section  III.2,  the  two  unit  power 
plane  wave  arrivals  impinging  on  the  array  have  incidence  angles  of  -  10°  and  -  20° .  The 

array  is  10  —  wavelength  long  and  has  32  sensors,  —  =  3.  Five  different  subarray  lengths 

are  selected:  28,  24,  20,  16,  14  so  that  the  number  of  averages  in  the  smoothing  increases 
from  5  to  19. 

As  the  number  of  average  increases  or  subarray  length  decreases,  the  covariance 
matrix  becomes  more  and  more  Toeplitz,  as  can  be  seen  on  Figures  in.6  to  HI.  10.  Using 

Equation  (3.4.8)  with  —  =  3,  60  =  -  10,  and  0!  =-20,  the  number  of  averages  necessary  to 

decorrelate  the  two  arrivals  is  K  ~  18  which  corresponds  to  a  subarray  length  of  15.  One 
checks  that  for  subarray  lengths  equal  to  14  and  16,  the  spatially  smoothed  matrix  is 
almost  Toeplitz,  with  constant  phase  and  contant  magnitude  on  the  upper  diagonals. 

III.4.3.2  Effects  on  the  Output  of  the  MVDR  Beamformer 

Theoretical  results  [Shan  1985]  provide  the  criteria  for  the  detection  of  correlated 
arrivals  using  spatial  smoothing  and  high  resolution  processors.  When  a  good  estimate  of 
the  power  is  required,  adaptive  beamforming  must  be  used,  such  as  the  MVDR  beam- 
former.  It  thus  is  of  interest  to  study  the  amount  of  loss  due  to  correlation  at  the  output  of 


35 


the  MVDR  beamformer  after  spatial  smoothing.  The  simulation  described  hereafter  com¬ 
pletes  the  already  known  qualitative  information  about  the  general  behaviour  of  the 
angular  spectrum  estimate  [Reddy  1987]  by  providing,  in  a  particular  case,  quantitative 
information  on  the  loss  due  to  correlation. 

As  before,  a  lOy  wavelength  array  with  32  sensors  receives  a  pair  of  correlated 

(p=  1)  unit  power  plane  wave  arrivals  in  the  set  with  incidence  angles  {0°,  5°,  10°,  15°, 
20° ,  25° ,  30° ,  35° ,  40° } .  A  background  sensor  noise  of  -  20  dB  is  also  included.  The  possi¬ 
ble  angular  spacings  are  multiple  of  5°  which  is  close  to  the  Rayleigh  resolution  of  the 
array.  For  each  pair  of  arrival  combinations,  spatial  smoothing  followed  by  the  MVDR 
beamforming  is  performed  and  the  loss  corresponding  to  the  arrivals  measured.  The  loss 
is  the  difference  in  dB  between  the  true  power  (in  this  case  0  dB)  and  the  power  of  the 
arrivals  at  the  output  of  the  beamformer. 

Figure  III.  11  shows  the  influence  of  the  arrival  pair  angular  spacing  on  the  loss 
due  to  correlation.  Panel  A  indicates  that  without  any  smoothing  the  MVDR  beamformer 
suffers  almost  complete  signal  cancellation.  The  other  panels  correspond  to  different 
subsegment  lengths  from  30  to  14,  the  number  of  averages  varying  from  3  to  19.  Each 
panel  indicates  the  amount  of  loss  due  to  correlation  that  can  be  expected  from  a  given 
number  of  averages  and  a  given  pair  angular  seperation.  Each  angular  seperation 
corresponds  to  a  number  of  combinaison  of  arrivals,  e.g.  there  is  a  10°  seperation  for  the 
following  pairs:  (0°  ,  10°),  (5°  ,  15°),  (10°  ,20°),  (15°  ,25°),  (20°  ,30°),  (25°  ,35°),  (30°  ,40°). 
A  small  subarray  with  14  sensors  allows  perfect  decorrelation  of  arrivals  as  close  as  10" 
and  limits  the  correlation  loss  to  a  few  dB  for  5°  spacing.  A  subarray  of  24  elements  does 
a  good  job  at  decorrelating  arrivals  seperated  by  at  least  10° . 

Figure  ID.  12  shows  the  same  results  of  loss  due  to  correlation  as  a  function  of  the 
number  of  averages.  For  a  given  scenario,  it  allows  one  to  study  the  impact  of  the  number 
of  averages  on  the  final  angular  spectrum  and  to  select  a  subarray  length  for  each  case.  It 


36 


indicates  that  considerable  smoothing  must  be  performed  when  the  two  arrivals  are  close 
together.  A  5°  seperation  requires  at  least  16  averages  for  the  MVDR  beamformer  to 
yield  a  minimum  amount  of  loss. 

In  the  case  of  a  10°  seperated  pair.  Figure  III.  13  studies  the  impact  under  several 
smoothing  conditions  of  the  pair  location  with  respect  to  the  broadside  direction.  The 
effect  of  moving  an  arrival  pair  from  broadside  (0°  ,10°),  toward  endfire  (30°  ,40°)  is 
small  since  the  variations  of  loss  are  at  most  of  the  order  of  2  dB,  as  the  pair  direction 
arrives  at  steeper  and  steeper  angles. 
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Figure  m.6:  Magnitude  and  phase  of  the  upper  diagonals  of  the  28  by  28  spatially 
smoothed  covariance  matrix.  Diagonals  are  identified  by  looking  at  their  length.  Matrix 
index  corresponds  to  element  index  along  the  diagonal.The  32  sensor  line  array  (25  m 
interelement  spacing)  receives  two  correlated  arrivals  (p  =  l)  at  -  10°  and  -  20"  incidence 
angles  at  20  Hz.  The  sound  speed  is  1498  m/s. 
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Figure  m.7:  Magnitude  and  phase  of  the  upper  diagonals  of  the  24  by  24  spatially 
smoothed  covariance  matrix.  Diagonals  are  identified  by  looking  at  their  length.  Matrix 
index  corresponds  to  element  index  along  the  diagonal.The  32  sensor  line  array  (25  m 
interelement  spacing)  receives  two  correlated  arrivals  (p  =  l)  at  -  I0a  and  -  20°  incidence 
angles  at  20  Hz.  The  sound  speed  is  1498  m/s. 
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Figure  m.8:  Magnitude  and  phase  of  the  upper  diagonals  of  the  20  by  20  spatially 
smoothed  covariance  matrix.  Diagonals  are  identified  by  looking  at  their  length.  Matrix 
index  corresponds  to  element  index  along  the  diagonal.The  32  sensor  line  array  (25  m 
interelement  spacing)  receives  two  correlated  arrivals  (p=  l)  at  -  10°  and  -  20s  incidence 
angles  at  20  Hz.  The  sound  speed  is  1498  m/s. 
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Figure  m.9:  Magnitude  and  phase  of  the  upper  diagonals  of  the  16  by  16  spatially 
smoothed  covariance  matrix.  Diagonals  are  identified  by  looking  at  their  length.  Matrix 
index  corresponds  to  element  index  along  the  diagonal.The  32  sensor  line  array  (25  m 
interelement  spacing)  receives  two  correlated  arrivals  (p  =  l)  at  - 10°  and  -  20°  incidence 
angles  at  20  Hz.  The  sound  speed  is  1498  m/s. 
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Figure  HI.  10:  Magnitude  and  phase  of  the  upper  diagonals  of  the  14  by  14  spatially 
smoothed  covariance  matrix.  Diagonals  are  identified  by  looking  at  their  length.  Matrix 
index  corresponds  to  element  index  along  the  diagonaLThe  32  sensor  line  array  (25  m 
interelement  spacing)  receives  two  correlated  arrivals  (p  =  I)  at  -  10"  and  -  20"  incidence 
angles  at  20  Hz.  The  sound  speed  is  1498  m/s. 
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Figure  m.ll:  Loss  at  the  output  of  the  MVDR  beamformer  after  spatial  smoothing  for  a 
correlated  arrival  pair.  The  pair  angular  spacing  varies  from  5C  to  40" .  A  given  spacing 
corresponds  to  several  combinations  of  arrivals  at  {0°,  5°,  10°,  ...  ,  40°}.  Panel  A 
corresponds  to  no  smoothing,  Panels  B,  C,  D,  E,  F,  G,  H,  I,  J  correspond  to  smoothing 
subsegment  lengths  {30, 28, 26, 24, 22, 20, 18, 16, 14). 


Figure  HI.  12:  Loss  at  the  output  of  the  MVDR  beamfonner  as  a  function  of  the  number 
of  averages  performed  by  spatial  smoothing.  Panels  A,  B,  C,  D,  E,  F,  G,  H  correspond  to 
correlated  (p=  l)  arrival  pairs  with  angular  separations  {5°,  10®,  15",  20",  25",  30",  35", 
40"). 
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Figure  HI.  13:  Impact  on  the  loss  at  the  output  of  the  MVDR  beamformer  due  the  loca¬ 
tion  from  broadside  of  a  10®  separated  pair.  Panel  A  corresponds  to  no  smoothing  and 
Panels  B,  C,  D,  E,  F,  G,  H,  I,  J  correspond  to  smoothing  subsegment  lengths  {30,  28,  26, 
24,22,20, 18, 16, 14). 
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m.5  Modified  Spatial  Smoothing 

HI.5.1  Description  and  Interpretation 

A  second  version  of  spatial  smoothing  has  been  proposed  [Williamsl988]  and 
uses  an  estimate  of  the  covariance  matrix  which  is  the  modified  covariance  matrix  of  the 
combined  forward/backward  linear  prediction  algorithm  for  autoregressive  spectral  esti¬ 
mation  [Burg 1967].  The  modified  covariance  method  is  known  to  perform  better  in  time 
series  autoregressive  spectral  estimation  in  producing  less  spurious  peaks  [Marplel987], 
The  modified  covariance  matrix  also  is  used  in  the  maximum  entropy  wavenumber  pro¬ 
cessing  for  linear  arrays  [Barnard  1982]  and  in  the  Kumaresan-Tufts  high  resolution 
method  of  direction  of  arrival  estimation  [Kumaresanl983]. 

Im  modified  spatial  smoothing,  the  covariance  matrix  of  the  array  data  xf  and  of 
the  reversed  data  xb  is  estimated.  Defining 

x/=[*o*i  *«]  (3.5.1) 

xj  =  [xH  *„_]  •  ■  •  x0]  (3.5.2) 

then,  the  modified  covariance  matrix  is 

Xf  xF+xb  x*j  (3.5.3) 

This  is  made  possible  by  the  special  geometry  of  equispaced  line  arrays.  As  in  the  origi¬ 
nal  spatial  smoothing,  the  modified  spatial  smoothing  does  the  average  of  the  modified 
covariance  matrices  estimated  on  overlapping  subarrays.  Modified  spatial  smoothing  can 
be  formulated  as  [Williamsl988] 

M -t 

=  v,  £  fcf  X;  Xf1  +  J  (X,  X"  )•  jl  (3.5.4) 

M-j  +  li=o  l  J 

where  *  is  the  complex  conjugate  operation  and  J  is  the  reflection  matrix  or  reverse 
matrix  given  by 
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m.5.2  Performance  of  the  Modified  Spatial  Smoothing 

It  is  shown  in  [Williams  198  8]  that  if  p  correlated  arrivals  impinge  on  the  array,  it 
is  sufficient  to  average  ^  subarrays  of  p  sensors  to  recover  the  rank  of  the  signal  covari¬ 
ance  matrix,  that  is  to  have  det  S  *  0.  This  result  holds  under  some  specific  constraints 
detailed  in  [Williamsl988].  With  modified  spatial  smoothing,  a  M  sensor  array  is  able  to 

detect  at  most  —  correlated  arrivals  [Williamsl988].  This  result  is  of  great  interest 

because  it  indicates  that  the  number  of  averages  necessary  to  allow  the  determination  of 
the  direction  of  arrivals  with  high  resolution  methods  is  half  what  the  original  spatial 
smoothing  requires  and  allows  the  so-called  increase  in  efficient  array  aperture. 

As  with  spatial  smoothing,  an  estimate  of  the  power  for  each  arrival  often  is  of 
interest.  Then,  one  can  ask  what  decorrelation  rate  modified  spatial  smoothing  achieves 
compared  to  the  original  spatial  smoothing. 

The  study  outlined  in  Section  III.4.3  and  detailed  in  [Reddy  1987]  is  repeated 
here  for  the  modified  spatial  smoothing.  One  considers  two  correlated  arrivals  impinging 
on  a  Af  sensor  array  referenced  to  its  center  as  in  [Takaol987].  The  full  aperture  plane 
wave  steering  vectors  are  given  by 


Af  = 


,  .2kA  M+ 1 .  .  „ 
expO-r“(* - JT-  )sm0i 


(3.5.6) 


The  array  covariance  matrix  can  be  seperated  into  a  part  that  corresponds  to  the  original 
spatial  smoothing  and  a  part  that  corresponds  to  the  modified  spatial  smoothing.  The  spa¬ 
tially  smoothed  part  of  the  signal  covariance  matrix  is  given  by 
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where 


Si  = 


M- 1 
2 
£ 


2(M-s+\)  .  _  m- l 
2 


D*  S  (  D*  )w 


(3.5.7) 


D  = 


e  x  0 


;2*Ycn& 


(3.5.8) 


0  c 

and  S  is  given  by  Equation  (3.3.8).  The  spatially  smoothed  array  covariance  matrix  is 
given  by 

R,  =  B  S,  Bh  (3.5.9) 

where  B  is  the  stacked  direction  matrix  over  a  subarray.  Carrying  out  the  algebra. 


S,(l, 2)  = 


1 


CT]  <J2  p 


sin^~-(sin  8q  -  sinGO 


sin-^-(.rin  60-sin9,) 
A, 


2(M-s+l) 

The  second  part  of  the  array  covariance  matrix  is  given  by 

R2  =  JRi‘J 


(3.5.10) 


(3.5.11) 


or 


R2  =  JB*S1’Bt  J  (3.5.12) 

Since  it  is  shown  in  [Williams  1988]  that  JB*  =B,  then  R2=BS*Bw  and  the  equivalent 
modified  pan  of  the  signal  covariance  matrix  is  given  by  §2  =  S *.  Then  the  (1,2)  term  of  S 
is 


S(1 , 2)  = 


g]CT2  Ip  I 
M-s+1 


sinj i~K(sinQ0  -  sinGj) 

- - - c.os4> 

sinx— (sm0o  -  sinQi) 


(3.5.13) 


where  <(>  and  Ip  I  are  given  by  the  correlation  between  the  two  arrivals  p=  Ipl  eJ*. 

This  result  shows  a  peculiarity  of  the  modified  spatial  smoothing:  when  <j>,  the 
electrical  phase  between  the  two  arrivals  at  the  center  of  the  array  (the  reference)  is  y, 
perfect  decorrelation  is  achieved.  The  array  covariance  matrix  is  Toeplitz  and  there  is  no 
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loss  due  to  correlation  at  the  output  of  the  MVDR  beamformer,  even  if  s  =  M .  When  the 
electrical  phase  is  zero,  the  decorrelation  rate  is  identical  to  that  of  spatial  smoothing. 
These  results  that  strongly  depend  on  the  electrical  phase  4>  are  similar  to  what  was 
observed  in  [Gabriell986]  for  adaptive  antennas  and  in  [Whitel979]  for  the  minimum 
entropy  method  operating  on  fully  correlated  arrivals  in  quadrature. 

The  chances  of  having  two  signals  in  quadrature  at  the  center  of  the  array  are 
slim.  Thus,  all  one  can  state  is  that  the  modified  spatial  smoothing  works  as  well  or  better 
than  the  original  spatial  smoothing  in  decorrelating  a  pair  of  correlated  arrivals. 

If  the  reference  is  taken  at  another  spatial  position  along  the  array,  the  depen¬ 
dence  of  the  modified  spatial  smoothing  performance  on  the  electrical  phase  <t>  still  exists. 
The  precise  phase  relationship  between  the  two  signals  for  which  decorrelation  is 

achieved  is  not  y  and  depends  on  the  reference  point,  the  interelement  spacing,  the 

wavelength,  and  the  arrival  angles. 

m.5.3  Simulation  Results 

To  illustrate  the  theoretical  developments  of  the  previous  sections,  the  effects  of 
modified  spatial  smoothing  with  two  signals  with  a  varying  electrical  phase  at  the  center 

of  a  32  element  array  now  is  studied.  As  before,  the  array  is  10  y  wavelength  long  and 

receives  unit  power  plane  wave  arrivals  at  -5°  and  0°  incidence  angle.  The  subarray 
length  is  32  so  that  only  two  covariance  matrix  averages  are  performed.  The  maximum 
loss  due  to  correlation  for  the  two  arrivals  at  the  output  of  the  MVDR  beamformer  is 
plotted  on  Figure  III.  14.  When  the  two  signals  are  in  quadrature,  there  is  no  loss  since  the 
modified  covariance  matrix  is  Toeplitz.  The  loss  due  to  correlation  increases  as  the 
electrical  phase  decreases,  and  the  loss  reduces  to  that  given  by  the  original  spatial 
smoothing,  on  the  order  of  50  dB.  Then,  the  behaviour  of  the  modified  spatial  smoothing 
is  identical  to  that  of  the  original  spatial  smoothing. 


The  angular  spectra  given  by  the  MVDR  beamformer  for  an  electrical  phase  of  0 
and  -j  radians  are  plotted  on  Figure  III.  15.  Perfect  resolution  is  achieved  when  the  two 
signals  are  in  phase  quadrature,  while  the  method  fails  when  the  two  signals  are  in  phase. 
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Figure  m.14:  Variations  of  the  loss  at  the  output  of  the  MVDR  beamfonner  after 
modified  spatial  smoothing  as  a  function  of  the  electrical  phase  between  two  correlated 
arrivals  (I  p  l  =  l)  with  -5°  and  0®  incidence  angles.  Phase  and  arrival  angles  are  relative 
to  the  aperture  center.  The  array  has  32  sensors  (25  m  spacing)  and  operates  at  20  Hz. 
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Figure  HI.  15:  MVDR  Angular  Spectra  using  the  full  array  modified  covariance  matrix  of 
two  correlated  arrivals  with  -5°  and  0°  incidence  angles.  In  Panel  A,  the  two  arrivals  are 
in  phase  (p=  1)  while  in  Panel  B  the  two  arrivals  are  in  quadrature  (p  =  exp;n/2)  with 
respect  to  the  array  center. 
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ni.6  Conclusions 

This  chapter  illustrates  the  effects  of  correlation  on  high  and  super  resolution 
techniques  and  explains  how  the  spatial  smoothing  preprocessing  techniques  work.  The 
smoothing  techniques  attempt  to  reduce  the  modulations  of  the  diagonals  of  the  covari¬ 
ance  matrix  which  should  be  Toeplitz  according  to  the  signal  model.  The  various 
interpretations  of  smoothing  were  summarized  and  their  performances  outlined. 

Guidelines  were  given  to  use  smoothing  with  the  MUSIC  algorithm  and  the 
MVDR  beamformer.  Spatial  smoothing  requires  d  averages  of  d  sensor  subarrays  for  the 
MUSIC  DOA  estimator  to  detect  d  correlated  signals,  while  the  modified  spatial  smooth¬ 
ing  generally  requires  averages.  The  different  nature  of  the  problem  that  the  MVDR 

beamformer  attempts  to  solve  was  pointed  out.  In  addition  to  the  direction  of  arrivals,  it 
estimates  the  power.  Signal  cancellation  occurs  as  soon  as  there  is  some  correlation 
between  arrivals.  The  difficulty  of  decorrelating  two  closely  spaced  arrivals  clearly 
appears  in  the  simulations  presented.  Some  substantial  loss  due  to  correlation  may  still 
occur  even  after  heavy  spatial  smoothing. 

Differences  between  the  original  spatial  smoothing  and  the  modified  spatial 
smoothing  were  outlined.  As  observed  in  the  simulations,  the  variations  of  the  decorrela¬ 
tion  rate  achieved  by  modified  spatial  smoothing  as  a  function  of  the  electrical  phase 
between  two  correlated  arrivals,  was  explained  analytically. 

Smoothing  reduces  the  effective  aperture  of  the  array,  and  therefore  the  resolu¬ 
tion.  In  addition,  it  increases  in  a  significant  way  the  computational  burden.  Nevertheless, 
it  is  a  necessary  preprocessing  step  before  using  the  high  or  super  resolution  techniques. 
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IV  Processing  Curved  Wavefronts 


The  effects  of  wavefront  curvature  are  studied  in  the  context  of  vertical  arrays 
using  a  geometric  ray  theory  approach.  Based  on  Snell’s  law,  a  way  to  generate  curved 
wavefront  is  derived  analytically.  The  importance  of  wavefront  curvature  first  is 
evaluated  by  comparing  the  phase  between  curved  wavefronts  and  plane  wavefronts, 
both  corresponding  to  the  same  arrival  angles  for  a  realistic  scenario.  The  behaviour  of 
the  conventional,  MVDR,  and  MUSIC  beamformers  using  plane  wavefront  replica  vec¬ 
tors  then  is  studied  in  a  curved  wavefront  environment.  The  largest  mismatch  loss  occurs 
for  near  horizontal  arrivals  and  the  aperture  length  is  not  critical  for  this  rather  general 
result.  Curved  wavefront  beamforming  then  is  proposed  and  the  behaviour  of  the  proces¬ 
sors  is  described  in  terms  of  bearing  responses  and  beam  patterns.  The  distortion  caused 
by  the  smoothing  transformation  is  studied  in  the  context  of  high  resolution  methods.  A 
method  using  the  eigenvectors  of  the  smoothed  covariance  matrix  for  each  replica  vector 
allows  the  loss  due  to  mismatch  to  be  avoided.  Finally,  the  performance  of  the  MVDR 
beamformer  after  spatial  smoothing  under  correlated  and  curved  arrivals  is  discussed. 
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IV.l  Introduction 

The  previous  chapter  began  a  realistic  description  of  the  oceanic  environment. 
The  importance  of  the  ocean  surface  and  bottom  was  noted  and  the  effects  of  multipath 
propagation  on  the  high  resolution  beamforming  methods  were  studied. 

Another  important  aspect  of  the  underwater  medium  is  its  highly  refractive  char¬ 
acter  which  is  variable  in  the  vertical  as  well  as  the  horizontal.  Propagation  over  long 
range  is  controlled  by  the  sound  velocity  profile  and  its  variations  in  range  [Mur- 
phyl987].  The  vertical  directionality  of  the  wavefield  is  a  natural  way  to  study  the  com¬ 
plex  propagation  mechanisms  due  to  the  environment,  and  vertical  line  arrays  are  major 
experimental  tools  used  to  measure  ambient  noise  as  well  as  long  range  propagation  sig¬ 
nals.  At  first  with  lengths  of  only  a  few  hundred  meters  [Anderson  1974,  Andersonl979, 
Kewleyl984,  Dossol987,  Senl988],  array  apertures  have  been  increased  in  recent  years 
to  provide  a  better  coverage  of  the  water  column  and  enable  work  at  lower  frequency 
[Sotirinl988].  In  this  framework,  it  is  of  interest  to  study  the  wavefront  curvature  due  to 
refraction. 

In  ray  theory,  the  variation  of  sound  speed  with  depth  causes  ray  bending  and 
can  be  summarized  by  the  Snell’s  law  [Brekhovskikhl982] : 

— =  constant  (4.1.1) 

cos0 

where  c  is  the  local  sound  speed  and  9  the  ray  angle  with  respect  to  the  horizontal.  Ray 
bending  can  be  large.  A  ray  with  angle  with  respect  to  the  horizontal  of  10°  at  great  depth 
with  sound  velocity  of  1500  m/s  has  an  angle  with  respect  to  the  horizontal  of  13.8°  shal¬ 
lower  in  the  water  column  where  the  sound  velocity  is  1480  m/s.  This  is  depicted  on  Fig¬ 
ure  HI.  1 

Since  ray  bending  is  large,  wavefront  curvature  is  likely  to  have  a  large  impact 
on  the  outputs  of  the  beamforming  structures.  Processing  a  field  with  varying  arrival 
angle  across  the  array  aperture  is  similar  to  the  case  of  processing  nonstationary  random 
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processes  in  time  series  analysis.  The  selection  of  the  plane  wave  array  manifold  that 
corresponds  to  a  medium  with  a  constant  sound  velocity,  may  not  be  appropriate,  espe¬ 
cially  for  a  very  long  vertical  line  array. 

Although  the  effects  of  index  of  refraction  variations  may  be  negligible  for  the 
conventional  processor  due  to  its  lack  of  sensitivity  or  its  robustness,  there  are  some  indi¬ 
cations  that  they  can  be  large  for  high  resolution  processors  [Seglisonl970,  Mac- 
Donoughl971]. 


IV .2  Generating  Curved  Wavefronts 


IV.2.1  Derivation  of  the  Phase  Relationships  in  the  Vertical 

In  order  to  understand  the  impact  of  curved  wavefronts  on  the  processing  struc¬ 
tures  described  in  the  previous  chapters,  a  model  of  the  wavefront  curvature  needs  to  be 
derived.  The  simple  geometric  approach  of  ray  theory  is  selected.  Since  the  sound  velo¬ 
city  continuously  varies  with  depth,  the  phase  variations  from  a  point  Ma  in  space  to  a 
point  M  j  can  be  expressed  by  the  integral  relationship 

*=/*</?  (4.2.1) 

M. 

where  is  the  phase  of  the  propagating  wavefront,  k  the  wavenumber  vector  varying 
over  space,  and  position  in  the  two  dimensional  space  sampled  by  the  vertical  array  is 
defined  by  the  coordinate  system  (r,z),  with  r,=  rel+ze ?.  (c7a)  spans  the  space  as 
defined  by  Figure  IV.2.  Fean  be  expressed  as 


*,  ,  2 nf 

( r,z  )=  7  - 

c(r,z) 


cos9(r,z)  e?  +  sin9(r,z)  e? 


(4.2.2) 


where  (r,  z)  are  the  coordinates  of  the  point  in  space  where  the  observation  is  made,  /  is 
the  acoustic  frequency,  c(r  ,z)  the  sound  velocity  at  the  point  (r  ,z)  and  9(r  ,z)  the  ray 
angle  with  respect  to  the  horizontal  at  (r ,  z). 
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The  receiving  array  has  its  sensors  numbered  from  top  to  bottom  and  negative 
angles  of  arrival  corresponds  to  downgoing  sound  or  uplooking  beams. 

In  the  context  of  vertical  line  arrays,  one  is  interested  in  the  phase  variations  due 
to  the  sound  speed  variations  in  the  vertical  so  that 

dr=dzel  (4.2.3) 

and,  assuming  a  stratified  ocean  (range  independant).  Equation  (4.2.1)  becomes 

4  =  2nf  ]  (4.2.4) 

i  c(*> 

Using  the  Snell’s  law,  one  gets 

<|>  =  2n/  J  5gw(e(z)}[  1  —  cos ^8(2 0)  }  dz  (4.2.5) 

i.  c(z)  l  c(Zo)J 

c(z)  2 

where  Sgn  is  the  signum  function,  and  cos2©^)  ;  ;  <  1  is  implicitely  assumed  for  all 

C(zo)J 

2.  If  the  medium  is  isovelocity,  i.e.  for  all  z  c(z)  =  c(z0 ),  Equation  (4.2.5)  reduces  to  the 
plane  wavefront  result : 


2nf  sin©(z0) 


<J\  -*o) 


(4.2.6) 


Equation  (4.2.5)  provides  the  variation  of  phase  in  the  vertical  and  is  identical  to 
a  ray  representation  of  the  field.  Following  [Ahluwalial977],  the  pressure  field  from  ray 
theory  is  given  by 


p(r  ,z)  =  expO^oS (r  ,z))  L  (jk0)m  Am(r  ,z)  (4.2.7) 

m«0 

The  summation  term  corresponds  to  the  amplitude  and  k0S(r  ,z)  gives  the  phase  varia¬ 
tions  in  space.  k0  is  the  reference  wavenumber  and  S(r  ,z)  often  is  being  called  the 
eikonal  [Brekhovskikhl982].  The  amplitude  and  the  eikonal  are  solutions  of  the 
Helmholtz  equation.  For  horizontally  stratified  ocean,  the  S(r  ,z)  can  be  expressed  as 
[Ahluwalial977,  pp  82-83] 
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i  <> 


v  r  t  * 

S(r,z)  =  ar  +  J  n2-a2J  dz 


(4.2.8) 


7  r  "i  % 
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dz 


(4.2.9) 


where  z0  is  the  reference  depth  taken  at  the  sound  axis,  n  is  the  refraction  index  given  by 
c(z0) 


n  = 


c(z) 


<  l,  and  a  is  a  constant  equal  to  n(z0)  sina  where  a  is  the  take-off  angle  in  spher¬ 


ical  polar  coordinate  (that  is  a=  y-0).  Then  Equation  (4.2.9)  is  equivalent  to  Equation 


(4.2.5). 


IV.2.2  Partial  Insonification 

The  ability  to  generate  a  curved  wavefront  array  manifold  will  allow  us  first  to 
simulate  a  realistic  field  consisting  of  curved  wavefronts,  and  second  will  constitute  a  set 
of  replica  vectors  for  use  with  the  high  resolution  structures. 

The  final  estimate  of  the  angular  spectrum  greatly  depends  on  the  selected  refer¬ 
ence  point,  that  is  the  depth  to  which  the  angular  spectrum  corresponds.  In  the  following, 
the  reference  is  chosen  at  the  center  of  the  array  aperture.  The  results  of  Section  IV.2.1 
are  used  to  compute  the  phase  entering  into  the  steering  vectors  stacked  in  the  array  man¬ 
ifold  matrix.  The  array  manifold  matrix  has  complex  columns  that  correspond  to  physical 
angles  from  -  90°  to  90°  and  lines  that  correspond  to  the  different  sensors  of  the  array. 
The  sensors  are  numbered  from  top  to  bottom.  The  i*  sensor  or  the  i*  line  in  the  matrix 
is  associated  with  a  sound  velocity  c, .  If  c,  >  c0,  where  c0  is  the  reference  sound  velocity 
(sound  velocity  at  the  center  of  the  array),  a  ray  with  angle  0O  with  respect  to  the  horizon¬ 
tal  at  the  center  may  not  propagate  to  the  sensor  depending  on  0O.  The  condition  for  phy¬ 
sical  existence  of  a  ray  at  the  i*  sensor  is 

0O  £  cos-1— 


(4.2.10) 
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If  Equation  (4.2.10)  is  not  verified,  the  ray  has  turned  over,  and  the  steering  vector  com¬ 
plex  exponential  is  replaced  by  zero.  Beyond  the  turning  point,  there  is  a  shadow  zone 
where  the  solution  of  the  wave  equation,  under  the  WKB  approximation,  is  exponentially 
decaying  [Boylesl984,  p  210].  Here,  the  model  adopted  assumes  that  there  is  no 
insonification  past  the  turning  point.  Figure  IV.  3  shows  the  array  manifold  stacked  steer¬ 
ing  matrix  for  a  particular  sound  velocity  profile  with  respect  to  a  receiving  line  array. 

IV.2.3  Phase  Difference  between  Curved  and  Plane  Wavefronts 

The  relative  importance  of  wavefront  curvature  needs  to  be  assessed.  This  can 
be  done  by  comparing  the  phase  of  a  curved  wavefront  to  that  of  a  plane  wavefront,  with 
the  reference  at  the  center  of  the  array  in  both  cases. 

A  simulation  is  performed  on  a  64  element  array  with  25  m  spacing  in  the  water 
column  with  a  sound  velocity  of  1481  m/s  at  the  top  of  the  array.  The  sound  velocity 
linearly  increases  across  the  array  to  1515  m/s  at  the  bottom  of  the  array.  The  frequency 
is  20  Hz.  The  phase  of  the  plane  wavefronts  is  computed  using  Equation  (4.2.6)  with  the 
sound  velocity  at  the  center  of  the  array,  1498  m/s,  as  reference  sound  speed. 

The  differences  of  phase  for  each  sensor  position  when  insonified  by  arrival 
angles  between  -  90°  and  90°  are  overlaid  on  Figure  IV. 4.  The  phase  change  across  the 
entire  array  for  the  endfire  arrival  is  on  the  order  of  66  radians.  One  therefore  notes  that 
the  difference  of  phase  varies  as  much  as  10  %  of  this  maximum  phase.  Figure  IV.4  also 
indicates  that  only  a  few  steering  vectors  yield  a  large  phase  difference  between  the 
curved  and  plane  cases.  The  low  angles  actually  are  the  only  ones  subject  to  these  large 
deviations,  and  this  effect  still  exists  when  the  array  aperture  is  reduced.  A  16  sensor 
aperture  from  sensor  number  24  to  40  will  experience  relatively  important  phase  varia¬ 
tions  at  low  angles. 


Phase  Difference  (rad) 


Figure  IV  .4:  Overlaid  plots  of  the  phase  difference  between  curved  and  plane  wavefronts 
corresponding  to  angles  with  respect  to  the  horizontal  between  -  90°  and  90° .  The  line 
array  has  64  sensors  (25  m  interelement  spacing)  in  a  constant  gradient  sound  velocitv. 
The  reference  sound  velocity  at  the  center  of  the  aperture  is  1498  m/s. 
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IV.3  Plane  Wave  Beamformers  in  a  Curved  Wavefront  Environment 

IV .3.1  Performance  of  the  Plane  Wave  Beamformers 

The  response  of  the  beamformers  to  a  curved  wavefront  field  is  of  interest  in 
order  to  assess  the  effects  of  the  curvature  on  their  performance. 

The  same  64  element  array  as  in  Section  IV.2.4  is  used  in  this  simulation.  The 
array  operates  at  20  Hz,  and  has  a  25  m  interelement  spacing.  The  sound  velocity  has  a 
constant  gradient  increasing  from  1481  m/s  at  the  top  of  the  array  to  1515  m/s  at  the  bot¬ 
tom.  The  data  at  the  input  of  the  plane  wave  beamformer  are  curved  wavefronts  gen¬ 
erated  as  discussed  in  the  previous  sections.  A  Kaiser-Bessel  window  with  an  a  parame¬ 
ter  of  1.5  weights  the  data  across  the  array  in  order  to  yield  a  38  dB  side  lobe  rejection. 
Bearing  responses  for  arrivals  at  -  40° ,  -  30° ,  -  20° ,  -  10° ,  -  5° ,  -  0.2° ,  0.2° ,  T ,  10° ,  20°  30° , 
40°  are  plotted  on  Figure  IV.5.  These  should  be  compared  to  the  bearing  responses  for 
plane  wavefront  data  which  are  plotted  on  Figure  IV.6. 

The  bearing  response  of  the  conventional  beamformer  for  curved  data  are  similar 
to  the  one  for  plane  data  in  case  of  high  angles  of  arrival.  Differences  clearly  appear  for 
low  angles  of  arrival  since  the  peak  indicating  the  curved  arrival  is  smeared.  The  estimate 
suffers  some  loss  for  arrival  angles  between  -  10°  and  10°,  and  some  large  angle  bias  for 
near  horizontal  arrivals  (e.g.  the  0.2®  arrival  yields  a  peak  at  5°  with  a  -  3  dB  loss).  Furth¬ 
ermore  the  side  lobe  level  is  raised  at  least  of  20  dB  in  the  curved  wavefront  case  with 
respect  to  the  plane  wave  case. 

As  indicated  in  Chapter  II,  the  MVDR  beamformer  has  a  high  sensitivity  to 
mismatch  compared  to  the  conventional  processor  [Seglisonl970,  McDonoughl971], 
The  plane  wave  MVDR  bearing  responses  are  plotted  on  Figure  IV.7  for  curved  data  and 
on  Figure  IV.8  for  plane  data.  With  curved  wavefront  data,  one  notes  that  the  MVDR 
processor  suffers  a  considerable  loss  due  to  mismatch.  In  addition,  the  bearing  response 
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for  arrival  angles  between  - 10°  and  10°  are  smeared  and  biased  in  angle. 

The  bearing  response  results  of  the  MUSIC  algorithm  are  similar  to  those  of  the 
MVDR  beamformer  results  as  shown  in  Figure  IV.9  (plane  plane  wave  MUSIC  process¬ 
ing  curved  wavefront  data)  and  Figure  IV.  10  (plane  wave  MUSIC  processing  plane 
wavefront  data). 

In  this  particular  case,  one  can  conclude  that  the  conventional  beamformer  some¬ 
what  suffers  from  an  imperfect  wavefront  modeling.  The  deleterious  effects  occur  where 
the  curvature  is  the  largest,  that  is  for  low  angle  arrivals.  The  effects  of  imperfect  wave- 
front  modeling  are  much  larger  for  the  MVDR  and  MUSIC  algorithms.  This  loss  due  to 
mismatch  cannot  be  afforded  and  curved  wavefront  replica  vectors  should  be  used. 

IV .3.2  Influence  of  Aperture  Length 

The  array  discussed  in  Section  IV. 3.1  is  1575  m  long  (or  approximately  21 
wavelengths  at  20  Hz)  and  the  question  of  the  dependence  of  curvature  effects  on  aper¬ 
ture  length  can  be  raised.  The  same  simulation  is  repeated  for  subarrays  extracted  from 
the  64  sensor  array.  These  subarrays  have  same  center  as  the  large  aperture  array.  The 
loss  due  to  mismatch  at  the  output  of  the  conventional  and  the  MVDR  beamformers  is 
plotted  on  Figure  IV.  1 1  for  the  different  subarray  lengths  and  for  arrival  angles  between 
-  40°  and  40° .  The  loss  due  to  mismatch  increases  as  the  subarray  length  increases  and  as 
the  arrival  becomes  close  to  the  horizontal.  While  the  curvature  has  only  a  mild  effect  on 
the  output  of  the  conventional  processor,  it  has  a  strong  effect  on  the  MVDR  results  for 
any  subarray  length  at  low  angles.  Figure  IV.  11  clearly  shows  that  curvature  cannot  be 
neglected  even  for  fairly  small  apertures  in  the  case  of  a  high  resolution  beamformer. 


Figure  IV.7:  Bearing  responses  of  the  MVDR  processor  using  plane  wave  replica  vectors 
in  a  curved  wavefront  field.  Panels  A,  B,  C,  D,  E,  F,  G,  H,  1,  j,  K,  L  correspond  respec¬ 
tively  to  arrivals  at  -  40° ,  -  30° ,  -  20° ,  -10° ,  -  5%  -  0.2° ,  0.2° ,  7" ,  10° ,  20° ,  30" ,  40" . 
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Figure  IV.9:  Bearing  responses  of  the  MUSIC  processor  using  plane  wave  replica  vec¬ 
tors  in  a  curved  wavefront  field.  Panels  A,  B,  C,  D,  E,  F,  G,  H,  I,  J,  K,  L  correspond 
respectively  to  arrivals  at  - 40° ,  -  30° ,  -  20",  -10“\  -  5° ,  -  0.2° ,  0.2° ,  T ,  10° ,  20" ,  30° ,  40" . 
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Figure  IV.  10:  Bearing  responses  of  the  MUSIC  processor  using  plane  wave  replica  vec¬ 
tors  in  a  plane  wavefront  field.  Panels  A,  B,  C,  D,  E,  F,  G,  H,  I,  J,  K,  L  correspond 
respectively  to  arrivals  at  - 40° ,  -  30° ,  -  2Xf , -\(f ,  -  5° ,  -  0.2° ,  0.2° ,  7",  10°,  20®,  30",  40°. 


Figure  IV.ll:  Loss  due  to  mismatch  for  the  conventional  and  MVDR  beamformers  using 
plane  wave  replica  vectors  under  curved  wavefront  data  for  several  array  lengths  (M  is 
the  number  of  sensors).  Panel  A:  conventional  processor.  Panel  B:  MVDR  processor. 
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IV.4  Curved  Wavefront  Beamforming 

IV .4.1  Conventional  Curved  Wavefront  Beamforming 

In  Chapter  II,  the  beamforming  operation  was  decomposed  into  two  components, 
a  processing  structure  and  an  array  manifold.  The  results  of  the  previous  section  indicate 
that  the  selection  of  the  plane  array  manifold  for  an  array  in  a  curved  wavefront  environ¬ 
ment  is  not  appropriate,  especially  for  high  resolution  methods  and  low  angle  of  arrivals. 
Curved  wavefront  steering  vectors  should  be  used  as  replica  vectors  [Murphy87]. 

The  Bartlett  beamformer  generally  is  used  with  a  tapering  window  in  order  to 
provide  good  sidelobe  rejection.  The  beampattems  of  the  untapered  Bartlett  beamformer 
steered  to  - 40°,  -  30°,  - 20° ,  -  10° ,  -  5° ,  - 0.2° ,  0.2° ,  7M0° ,  20° ,  30° ,  and  40°  are  plotted  on 
Figure  IV.  1 1.  The  first  sidelobes  is  around  -15  dB,  which  is  not  acceptable.  Windowing  is 
necessary  but  not  as  simple  as  in  the  plane  wave  case,  due  to  partial  insonification  of  the 
array  aperture.  The  window  length  has  to  be  adjusted  and  applied  to  the  insonified  part  of 
the  array  manifold  stacked  steering  matrix.  In  Figure  IV.  13  are  plotted  the  beam  patterns 
when  a  Kaiser-Bessel  window  with  an  a  parameter  of  1.5  has  been  properly  applied.  The 
sidelobe  rejection  is  equivalent  to  the  plane  wave  case  except  at  low  angles  where  there 
is  a  large  sidelobe  around  -  30  dB 

This  high  sidelobe  indicates  that  the  beair former,  steered  in  the  40°  direction, 
still  receives  energy  from  the  near  horizontal  attenuated  by  30  dB.  This  can  be  also 
observed  in  Figure  IV.  14  where  bearing  responses  for  the  same  angles  are  plotted.  The 
side  lobe  floor  is  raised  when  the  arrival  angle  approaches  the  horizontal.  This  sidelobe  is 
due  to  the  partial  insonification  of  the  aperture. 


Figure  IV.  12:  Beampatterns  of  the  conventional  beamformer  using  curved  wavefront 
replica  vectors  under  curved  wavefront  data.  Panels  A,  B,  C,  D,  E,  F,  G,  II,  I,  J,  K,  L 
correspond  to  steering  angles  {-40°,  -30°,  -20°,  -10°,  -5°,  -0.2°,  0.2®,  7®,  10®,  20®,  30®, 
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IV .4.2  High  Resolution  Curved  Wavefront  Beamforming 


IV .4.2.1  Effects  of  Curvature  on  the  Array  Covariance  Matrix 

To  understand  why  the  plane  wave,  high  resolution  beamformers  do  so  poorly 
under  curved  wavefront  data,  the  effects  of  curvature  on  the  array  covariance  matrix  are 
now  investigated.  Uncorrelated  plane  wavefront  arrivals  impinging  on  a  line  array  with 
equispaced  sensors  lead  to  a  Toeplitz  Hermitian  covariance  matrix.  When  the  arrivals 
have  curved  wavefronts,  the  corresponding  phase  variations  from  sensor  to  sensor  are  not 
constant  as  they  are  for  plane  waves.  If  A  is  a  model  signal  direction  matrix  expressed  as 

Ar  =  [e'*e'*  •■•«'*]  (4.5.1) 

its  contribution  to  the  array  model  covariance  matrix,  a  aw,  is  a  matrix  with  first  upper 
diagonal  given  by 


and  because  of  curvature,  one  has  for  i*j 


(4.5.2) 


*  Qj  ~  fy+i  (4.5.3) 

so  that  A  A"  is  not  Toeplitz  and  the  model  covariance  matrix  is  not  as  well. 

Figure  IV.  15  and  Figure  IV.  16  illustrate  the  non-Toeplitz  character  of  the  covari¬ 
ance  matrix  in  the  case  of  a  single  arrival  with  a  -  10"  angle  with  respect  to  the  horizon¬ 
tal,  and  an  uncorrelated  arrival  pair  with  -  20°  and  -  10"  arrival  angles.  The  scenario  is 
the  same  as  before  with  a  32  element  array  extracted  from  the  large  64  sensor  array  of 
Section  IV.2.4  and  centered  at  the  same  depth.  Curvature  makes  Figure  IV.  16  look 
significantly  different  from  Figure  III. 3,  its  plane  wave  counterpart. 
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Figure  IV.  15:  Upper  diagonals  of  the  32  sensor  array  covariance  matrix  of  a  unit  power 
curved  wavefront  arrival  with  -  10°  arrival  angle.  Matrix  index  corresponds  to  index  on 
the  diagonal.  Each  diagonal  can  be  identified  by  its  length. 
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Figure  IV.  16:  Upper  diagonals  of  the  32  sensor  array  covariance  matrix  of  two  unit 
power  uncorrelated  curved  wavefront  arrivals  with  -20°  and  - 10°  arrival  angles.  Each 
diagonal  can  be  identified  by  its  length. 
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IV .4.2.2  Curvature  and  Spatial  Smoothing 

Under  correlated  arrivals,  the  high  resolution  methods  fail  unless  spatial  or 
modified  spatial  smoothing  is  performed.  It  was  shown  in  Chapter  III  that  these  prepro¬ 
cessing  techniques  specifically  are  designed  for  equispaced  sensor  arrays  in  a  plane 
wavefront  field.  They  use  the  special  underlying  structure  of  the  covariance  matrix:  a 
Toeplitz  structure. 

In  a  curved  wavefront  environment,  the  direct  application  of  these  preprocessing 
schemes  creates  a  twofold  problem.  One  is  the  distortion  introduced  by  the  smoothing 
procedure,  the  underlying  structure  of  the  covariance  is  not  Toeplitz  and  depends  on  the 
arrivals  which  are  unknown  a  priori.  The  other  is  related  to  the  actual  decorrelation  rate 
achieved  by  the  smoothing  procedure,  thus  directly  influencing  the  performance  of  the 
high  resolution  beamformers. 

TV.4.2.3  Limiting  the  Distortion  Effects  due  to  Smoothing 

The  study  of  distortion  created  by  smoothing  the  original  covariance  matrix  is 
based  on  the  bearing  response  of  the  MVDR  beamformer  and  the  MUSIC  processor  for 
arrival  angles  at  {  -40°,  -30°,  -20°,  -10°,  -5°,  -0.2°,  0.2°,  7°,  10°,  20°,  30°,  40°  }.  The 
simulation  is  done  under  the  same  conditions  as  Section  IV. 3.1  with  a  smaller  aperture  of 
32  seniors  extracted  from  the  large  64  sensor  array  with  the  same  center. 

After  smoothing  with  a  subarray  length  of  s  sensors,  a  natural  choice  for  the  set 
of  replica  vectors  is  the  curved  wavefront  array  manifold  that  corresponds  to  a  s  sensor 
array  extracted  from  the  large  aperture  with  same  center.  The  bearing  responses  of  the 
MVDR  beamformer  after  spatial  smoothing  with  subarray  lengths  of  24  and  16  respec¬ 
tively,  are  plotted  on  Figure  IV.  17  and  Figure  IV.  18.  The  bearing  responses  of  the 
MUSIC  algorithm  for  smoothing  subsegment  length  of  24  and  16  respectively,  are  plot¬ 
ted  on  Figure  IV.  19  and  Figure  IV. 20.  One  notes  the  loss  due  to  mismatch  for  arrivals  at 
low  angle. 
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This  loss  is  due  to  mismatch  between  the  curved  wavefront  replica  vector  and 
the  equivalent  array  signal  vector  after  spatial  smoothing.  It  is  the  result  of  the  distortion 
due  to  smoothing.  For  each  bearing  characterized  by  a  signal  direction  vector  A  and  the 
corresponding  covariance  matrix  A  Aw,  this  mismatch  loss  can  be  minimized  by  choosing 
the  replica  vector  E  that  maximizes  the  power  at  the  output  of  the  conventional  processor 


E " 


T (A  \H)  ^  E 


(4.2.4) 


where  T  denotes  the  transformation  performed  by  smoothing  on  the  covariance  matrix. 
The  best  replica  vector  is  given  by  the  eigenvector  of  T(\AH)  that  has  the  largest  eigen¬ 
value.  This  replica  vector  will  be  referred  to  as  the  eigen  steering  vector.  Because  the  dis¬ 
tortion  effects  associated  with  the  smoothing  transformation  are  relatively  small,  the 
smoothed  matrix  is  of  rank  1  for  all  practical  purposes.  An  eigenvalue/eigenvector 
decomposition  of  the  smoothed  version  of  the  (32,  32)  covariance  matrix  for  both  subar¬ 
ray  lengths  was  performed  for  256  steering  vectors  spanning  the  grazing  angles  from 
-90°  to  90°.  The  eigenvalues  are  overlaid  on  Figure  IV.21  for  the  subsegment  lengths  of 
24  and  16  respectively.  There  is  a  single  large  eigenvalue  of  value  one,  the  others 
decrease  to  the  normalized  fraction  of  noise  used  to  stabilized  the  matrix. 


Figure  IV.22  and  Figure  IV.23  give  the  bearing  responses  for  the  MVDR  beam- 
former  after  spatial  smoothing  with  subsegment  length  of  24  and  16  using  the  eigen 
steering  vectors.  Figure  IV.24  and  Figure  IV.25  give  the  bearing  responses  for  the 
MUSIC  algorithm  after  spatial  smoothing  with  subsegment  lengths  equal  to  24  and  16, 
using  the  eigen  steering  vectors.  Loss  due  to  mismatch  is  reduced  to  negligible  values. 
Nevertheless,  the  use  of  the  eigen  steering  vector  cannot  reduce  the  -  25  dB  side  lobe 
level  at  near  horizontal  angle  in  the  case  of  the  MVDR  beamformer  after  smoothing  with 
a  subsegment  length  of  16. 
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IV.20:  Bearing  responses  of  the  MUSIC  beamformer  on  a  32  element  array  after 
smoothing  with  a  subsegment  length  of  16  with  curved  wavefront  replica  vectors, 
ta  are  curved.  Panels  A,  B,  C,  D,  E,  F,  G,  H,  I,  J,  K,  L  correspond  to  arrivals 
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Figure  IV.21:  Overlaid  eigenvalues  of  the  spatially  smoothed  covariance  matrix  for  all 
256  curved  wavefront  steering  vectors  corresponding  to  arrival  angles  between  -  90°  and 
90°.  Panel  A  correponds  to  a  subarray  length  of  24,  while  Panel  B  to  a  length  of  16. 


gure  IV .23:  Bearing  responses  of  the  MVDR  beamformer  after  spatial  smoothing  with 
i  sensor  subarrays  using  the  eigensteering  vectors.  Panels  A,  B,  C,  D,  E,  F,  G,  1 1, 1,  J,  K, 
correspond  to  arrivals  {-  40° ,  -  30" ,  -  20" ,  -10" ,  -  5" ,  -  0.2" ,  0.2" ,  7" ,  10" ,  20" ,  30" ,  40" } . 


Figure  IV.24:  Bearing  responses  of  the  MUSIC  processor  after  spatial  smoothing  with 
24  sensor  subarrays  using  the  eigensteering  vectors.  Panels  A,  B,  C,  D,  F,  F,  G,  II,  I,  J, 
K,  L  correspond  to  arrivals  (-40°,  -30",  -2(F,  -10°,  -5",  -0.2°,  0.2°,  7°,  10°,  20°,  30°, 


,ure  IV.25:  Bearing  responses  of  the  MUSIC  processor  after  spatial  smoothing  with 
sensor  subarrays  using  the  eigensteering  vectors.  Panels  A,  B,  C,  D,  E,  F,  G,  II,  I,  J, 
L  correspond  to  arrivals  -  40° ,  -  30° ,  -  20° ,  -10° ,  -  5" ,  -  0.2" ,  0.2 0 , 7° ,  10" ,  20° ,  30" ,  40"  } . 
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IV .4.2.4  Performance  of  the  MVDR  Beamformer  after  Smoothing 

In  a  plane  wave  environment,  Chapter  ID  showed  that  the  performance  of  the 
MVDR  beamformer  after  modified  spatial  smoothing  is  similar  to  the  performance  after 
spatial  smoothing.  The  use  of  the  modified  covariance  matrix  in  a  curved  Advefront 
environment  is  not  easily  justifiable  since  the  symetry  present  in  the  plane  wavefront 
problem  is  lost  in  a  curved  wavefront  environment  Therefore,  this  study  of  performance 
will  be  only  concerned  with  spatial  smoothing. 

In  this  simulation,  the  same  32  element  array  receives  a  pair  of  unit  power  fully 
correlated  (p  =  1)  curved  wavefront  arrivals  with  arrival  angles  in  the  set  {  0" ,  5° ,  10° ,  15° , 
20° ,  25" ,  30° ,  35" ,  40°  } .  A  background  sensor  noise  of  -  20  dB  is  also  included.  Spatial 
smoothing  followed  by  the  MVDR  beamformer  is  performed  using  the  eigen  steering 
vectors  that  correspond  to  a  particular  smoothing  scenario.  The  loss  suffered  by  the 
arrivals  is  measured  by  taking  the  difference  of  the  true  underlying  power  (0  dB)  and  the 
power  reported  by  the  MVDR  beamformer.  The  results  summarized  in  Figure  IV.26  and 
Figure  IV.27.  Figure  IV.26  shows  the  influence  of  the  arrival  pair  angular  spacing  on  the 
loss  due  to  coherence.  Each  angular  separation  corresponds  to  a  number  of  combinations 
of  arrivals  (e.g.  the  5"  pair  spacing  corresponds  to  the  pairs  (0"  ,  5"),  (5" ,  10"), ...).  Figure 
IV.27  gives  the  same  results  in  term  of  number  of  averages.  This  simulation  correponds 
to  the  plane  wave  results  summarized  in  Figure  III.  1 1  and  Figure  III.  12  of  Chapter  HI. 

Panel  A  of  Figure  IV.26  indicates  that  without  smoothing  the  MVDR  suffers 
almost  total  signal  cancellation  as  in  the  plane  wave  case.  The  other  panels  correspond 
to  different  subsegment  lengths  from  30  to  14  (the  number  of  averages  varies  from  3  to 
19).  The  beamformer  fails  to  resolve  the  arrival  pair  (0",  5")  for  any  subsegment  length. 
When  one  of  the  arrivals  is  near  horizontal,  the  power  estimate  suffers  greater  loss  than 
for  a  pair  with  same  angular  separation  but  with  higher  arrival  angle  (Figure  IV.27).  A 
near  horizontal  arrival  has  the  most  curvature  and  partially  insonifies  the  array  (because 
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the  sound  velocity  increases  with  the  sensors,  only  the  upper  half  of  the  array  is 
insonified).  One  can  note  that  with  the  exception  of  this  case,  the  results  are  very  similar 
to  those  for  plane  waves. 
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Figure  IV.26:  Loss  at  the  output  of  the  MVDR  beamformer  after  spatial  smoothing  for  a 
correlated  pair  with  angular  separation  varying  between  5°  and  40° .  The  data  are  curved 
and  the  replica  vectors  are  the  eigensteering  vectors.  Panel  A  corresponds  to  no  smooth¬ 
ing  and  Panels  B,  C,  D,  E,  F,  G,  H,  I,  J  correspond  to  smoothing  subsegment  lengths  { 30, 
28,26,24,22,20,18,16,14}. 
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Figure  IV.27:  Variation  of  the  loss  at  the  output  of  the  MVDR  beamformer  as  a  function 
of  the  number  of  averages  in  the  spatial  smoothing  for  a  correlated  pair  with  angular 
separation  varying  from  5®  to  40°.  The  data  are  curved  and  the  replica  vectors  are  the 
eigensteering  vectors.  Panel  A,  B,  C,  D,  E,  F,  G,  H  correspond  to  pair  separations  {5°, 
10" ,  15° ,  20° ,  25® ,  30® ,  35® ,  40®  } . 
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IV.5  Conclusions 

The  first  part  of  this  Chapter  was  devoted  to  deriving  an  expression  for  curved 
wavefront  steering  vectors.  The  simple  geometric  approach  of  ray  theory  provides  a  con¬ 
venient  way  to  generate  a  curved  wavefront  array  manifold.  The  implemented  integral 
equation  of  phase  reflects  the  variation  of  sound  speed  with  depth. 

The  importance  of  wavefront  curvature  has  been  assessed.  It  was  shown  that  cur¬ 
vature  is  large  for  near  horizontal  arrivals,  even  for  a  relatively  short  aperture.  On  the 
other  hand,  arrivals  with  steep  angles  with  respect  to  the  horizontal  are  quite  similar  to 
plane  wave  arrivals.  This  result  is  observed  on  the  bearing  responses  of  the  various 
beamformers  using  the  plane  wave  array  manifold  with  good  results  for  steep  angles  with 
respect  to  the  horizontal  and  poor  results  at  near  horizontal  arrival  angles.  One  also 
observes  the  great  sensitivity  to  mismatch  of  the  high  and  super  resolution  techniques 
with  losses  as  large  as  40  dB  when  curved  wavefronts  arrivals  impinge  on  the  array.  To 
avoid  loss  due  to  mismatch,  curved  wavefront  replica  vectors  should  be  used. 

The  direct  application  of  spatial  smoothing  to  decorrelate  correlated  curved 
wavefront  arrivals  causes  some  slight  distortion,  still  yielding  a  loss  due  to  mismatch 
even  when  using  the  curved  wavefront  array  manifold.  The  loss  due  to  mismatch  can  be 
avoided  by  taking  the  replica  vectors  that  correponds  to  the  steering  vectors  after  the  spa¬ 
tial  smoothing  transformation  of  the  covariance  matrix.  By  using  these  eigen  steering 
vectors,  loss  due  to  mismatch  is  made  negligible.  The  performance  of  the  MVDR  beam- 
former  under  correlated  curved  wavefront  arrivals  then  was  shown  to  be  very  similar  to 
the  plane  wave  case  (plane  wave  data  and  plane  wave  replica  vectors),  except  for  near 
horizontal  arrivals  where  performance  degradation  still  occurs. 
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